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Abstract. We use absolute trigonometric parallaxes from 
the Hipparcos Catalogue to determine individual dis- 
tances to members of the Hyades cluster, from which 
the 3-dimensional structure of the cluster can be derived. 
Inertially-referenced proper motions are used to rediscuss 
distance determinations based on convergent-point anal- 
yses. A combination of parallaxes and proper motions 
from Hipparcos, and radial velocities from ground-based 
observations, are used to determine the position and ve- 
locity components of candidate members with respect to 
the cluster centre, providing new information on cluster 
membership: 13 new candidate members within 20 pc of 
the cluster centre have been identified. Farther from the 
cluster centre there is a gradual merging between certain 
cluster members and field stars, both spatially and kine- 
matically. Within the cluster, the kinematical structure is 
fully consistent with parallel space motion of the compo- 
nent stars with an internal velocity dispersion of about 
0.3 km s _1 . The spatial structure and mass segregation 
are consistent with A-body simulation results, without 
the need to invoke expansion, contraction, rotation, or 
other significant perturbations of the cluster. The qual- 
ity of the individual distance determinations permits the 
cluster zero-age main sequence to be accurately modelled. 
The helium abundance for the cluster is determined to be 
Y = 0.26 ±0.02 which, combined with isochrone modelling 
including convective overshooting, yields a cluster age of 
625 ± 50 Myr. The distance to the observed centre of mass 
(a concept meaningful only in the restricted context of the 
cluster members contained in the Hipparcos Catalogue) 
is 46.34 ± 0.27 pc, corresponding to a distance modulus 



* Based on observations made with the ESA Hipparcos 
astrometry satellite. Table 2 is also available in 
electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u- 
strasbg.fr/Abstract.html 



m—M = 3.33±0.01 mag for the objects within 10 pc of the 
cluster centre (roughly corresponding to the tidal radius). 
This distance modulus is close to, but significantly better 
determined than, that derived from recent high-precision 
radial velocity studies, somewhat larger than that indi- 
cated by recent ground-based trigonometric parallax de- 
terminations, and smaller than those found from recent 
studies of the cluster convergent point. These discrepan- 
cies are investigated and explained. 
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1. Introduction 

The considerable importance of the Hyades cluster in 
studies of Galactic structure, in the understanding of the 
chemical evolution of the Galaxy, and in the determi- 
nation of the Population I distance scale, is well doc- 
umented in the literature. The nearest moderately rich 
cluster, with some 300 possible members, a total mass of 
some 300-400 M , and an age of around 600-800 Myr, 
it has an extension in the sky of about 20 degrees. Al- 
though uncertainty in the distances of individual members 
has limited the definition of the cluster's main sequence, 
and thereby its helium content and corresponding evo- 
lutionary sequence, it has nevertheless been used as the 
basic observational material for several fundamental re- 
lationships in astrophysics, including the location of the 
main sequence in the Hertzsprung-Russell diagram and 
the mass-luminosity relationship, as well as forming the 
basis for the determination of luminosities of supergiants, 
OB stars, and peculiar stars in clusters. Determinations 
of the distance to the cluster have provided the zero-point 
for distances within our Galaxy and, indirectly through 
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the Cepheids, one of the foundations on which the extra- 
galactic distance scale ultimately rests. 

At 40-50 pc, the Hyades cluster is somewhat beyond 
the distance where the parallaxes of individual stars are 
easily measured, or generally considered as fully reliable, 
from ground-based observations. Over almost a century, 
considerable effort using a wide variety of indirect meth- 
ods has therefore been brought to bear on the problem 
of establishing the distance to the cluster. Distance esti- 
mates have been based on a variety of geometrical mani- 
festations of a cluster of stars participating in a uniform 
space motion, while other estimates have been based on 
the average trigonometric parallax for a number of cluster 
stars, dynamical parallaxes for binaries, and photometric 
parallaxes using a variety of photometric systems. Nev- 
ertheless, the details of the HR and mass-luminosity dia- 
grams remain imprecisely established due to limitations in 
the accuracy of the parallaxes of the individual members, 
while the distance of the cluster is still open to debate: 
recent estimates of the distance modulus range from 3.16 
based on trigonometric parallaxes (Gatewood et al. 1992), 
3.40 based on convergent-point analyses using proper mo- 
tions from the FK5 and PPM Catalogues (Schwan 1991), 
and 3.42 based on recent Hubble Space Telescope FGS 
observations (van Altena et al. 1997a). 

The present work is inspired by the availability of the 
final results of the Hipparcos astrometry mission, which 
provide a radical improvement in astrometric data on 
all stars in the Hipparcos observing programme, includ- 
ing approximately 240 candidate Hyades members. The 
Hipparcos results offer the following principal improve- 
ments: (1) standard errors of the annual proper motions 
of typically 1 milliarcsec (mas) with respect to an inertial 
(extragalactic) reference frame; (2) absolute trigonomet- 
ric parallaxes with standard errors of order 1 mas; (3) sys- 
tematic errors of the astrometric parameters below around 
0.1 mas (or mas/yr); (4) parameter determination, or indi- 
cations, of double or multiple systems for component sep- 
arations larger than about 0.1 arcsec and Am Ss 3 mag; 
(5) precise photometry and detailed variability indicators 
based on the Hipparcos broad-band magnitude Hp; (6) 
homogeneous B — V and V — I colour indices. 

Literature on the Hyades distance determination is 
considerable: this paper is intended neither as a compre- 
hensive review nor a critical evaluation of the previous 
estimates in the light of the Hipparcos results. Neither 
does it aim to answer unambiguously the question 'what 
is the distance to the Hyades', a somewhat nebulous prob- 
lem given the resolution in radial distance provided by 
the Hipparcos parallaxes, and the sensitivity of the results 
to the precise qualification of distance: if the distance of 
the centre of mass is the objective, membership criteria, 
selection effects, and M/L relationships become critical. 
Rather, our objective is to reconcile previous distance es- 
timates based on the availability of reliable absolute tri- 
gonometric parallaxes, assign improved membership prob- 



abilities, and thereafter probe both the cluster dynamics 
and the assumptions on which previous distance deter- 
minations have rested. Finally, we will define the obser- 
vational main sequence based on a subset of objects for 
which membership is secure and observational data par- 
ticularly reliable, and compare this with theoretical deter- 
minations of the Hyades zero-age main sequence based on 
knowledge of the cluster's metallicity. 

In order to establish the complexities of the problem, 
Sect. 2 provides a summary of (some of) the distance 
determinations discussed in the literature to date, with 
particular reference to the agreement or disagreement be- 
tween the results of the various convergent-point analyses 
and distance estimates derived by other means. Sect. 3 
summarises the data, both from the Hipparcos Catalogue 
and from the published literature, used for the present 
study. Selection effects entering the list of candidates for 
this study are also discussed. 

The development of the paper then proceeds as fol- 
lows. In Sect. 4, we examine the improvement brought by 
the Hipparcos proper motions (and their connection to an 
inertial frame) which, as we shall demonstrate, permit a 
significant advance in the understanding of the systematic 
effects entering previous evaluations of the distance to the 
Hyades based on convergent-point analyses. New insights 
and the limitations of this approach applied to the Hyades 
are discussed. In Sect. 5 we use the Hipparcos absolute tri- 
gonometric parallaxes to determine a statistically signifi- 
cant distance estimate for each candidate member, even- 
tually permitting a provisional mean cluster distance to 
be defined. This is carried out in parallel with a combina- 
tion of the parallaxes and proper motions with published 
and unpublished radial velocities to determine the posi- 
tion and velocity components of candidate members with 
respect to a reference cluster centre. In Sect. 6 we discuss 
the Hipparcos parallaxes: first we combine the information 
coming from the Hipparcos parallaxes and proper motions, 
and demonstrate their mutual consistency. Then we exam- 
ine the differences between ground-based and Hipparcos 
parallaxes. Finally, we examine effects (in particular 'Lutz- 
Kelker' type corrections) which complicate the direct in- 
terpretation of the Hipparcos parallaxes. 

In Sects 7 and 8 we examine the spatial distribution 
and dynamics of the cluster, looking at the question of 
mass segregation, and comparing our present results with 
published A^-body simulations. We examine the velocity 
residuals of each member, which can be fully explained on 
the basis of the observational errors. We also examine the 
consistency of these results with previous estimates of the 
rotation, flattening, and internal velocity dispersion of the 
system. 

Finally, in Sect. 9, we restrict our list of Hyades can- 
didates to those showing no existing evidence for mul- 
tiplicity, and construct the resulting observational main- 
sequence, comparing it with new models for the zero-age 
main sequence. From this, the cluster helium abundance 
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is estimated, and this information is used to construct 
isochrones from which an estimate of the cluster age is 
determined. 



2. Distance determinations to date 

Under certain conditions which are at least reasonably well 
represented in the Hyades, the common proper motion 
of stars in a cluster can be used in the well-established, 
but nonetheless ingenious 'convergent-point' method of 
distance determination. Since all distance determinations 
employed for the Hyades have either been based on this 
method, or have been judged according to their agreement 
or otherwise with it, a short review of its theoretical basis 
and a comparison between its results and those obtained 
by other methods over the years is in order. Reviews of 
the various distance determination methods include those 
by van Altena (1974), Hanson (1980), and Turner et al. 
(1994). We will not discuss details of distance estimates 
based on photometric parallaxes (see, e.g., van Altena 
1974, Turner et al. 1994). A summary is given in Table 1. 

,/From simple dynamical arguments it can be shown 
that, for an open cluster of a few hundred stars within a 
volume of a few parsecs in radius, moving together under 
their mutual gravitation, the internal velocity dispersion 
is of the order of 1 km s _1 or less, and thus small com- 
pared with the typical linear velocity of the cluster as a 
whole relative to the Sun, v (the Hyades cluster has a 
space motion of approximately 45 km s^" 1 with respect to 
the Sun, presumably reflecting the velocity of the cloud 
in which the cluster formed). Provided that the cluster is 
sufficiently nearby to extend over an area of, say, several 
degrees, the parallel motions of the stars in space yield, 
on the celestial sphere, directions of proper motions that 
appear to converge on a unique point - the direction of 
the unit vector (v) is known as the convergent point. If b 
is the barycentric coordinate vector to a cluster member, 
and (b) its coordinate direction then, neglecting the in- 
ternal velocity dispersion, the radial velocity is p = (b)'v 
(the prime symbol associated with matrices and vectors 
will be used to denote transposition, so that x'y denotes 
the scalar product of the two vectors). With A denoting 
the angular distance between the star and the convergent 
point, and n the proper motion vector: 



(b)'(v) = cos A 
and: 

\n\ = 7T |v| sin X/A v 



(1) 



(2) 



where it is the parallax of the cluster member. A v = 
4.74047... km yr s _1 is the astronomical unit expressed 
in the appropriate form when tt and /x are expressed in 
mas and mas yr -1 respectively. For the Hyades, A ~ 33°, 
the radial velocity w rac j — 40 km s _1 in the cluster centre, 
|v| ~ 45 km s -1 , and ~ 100 mas yr -1 . 



Table 1. Distance determinations ordered by date (... in- 
dicates one or more co-authors). The distance modulus 
(D.M.) is taken or derived from the original reference, oc- 
casionally with some uncertainty as to a definitive value 
(and not necessarily in agreement with values referred to 
in subsequent compilations). 



D.M. 


Year 


Author 


Method 


2.75 


1939 


Smart 


Convergent point 


2.91 


1945 


Scares 


Convergent point 


3.03±0.06 


1952 


van Bueren 


Convergent point (GC, N30, NWZC) 


2.85 


1955 


Pearce 


Convergent point 


3.08 


1956 


Hcckmann... 


Convergent point 


3.04 


1965 


Wayman... 


Convergent point 


3.23±0.12 


1967 


Wallcrstein... 


Dynamical parallaxes 


3.14±0.19 


1967 


Eggen 


Trigonometric parallaxes (Yale) 


3.37 


1967 


Iben 


Stellar interiors 


3.08±0.07 


1969 


Sears... 


Stebbins photometric parallaxes 


3.10±0.06 


1969 


Eggen 


R — I photometric parallaxes 


3.25±0.20 


1969 


Heifer 


Wilson-Bappu 


3.09±0.06 


1970 


Upton 


pm gradient (FK4, N30, Yale, SAO) 


3.23±0.25 


1970 


Lutz 


Wilson-Bappu 


3.19±0.06 


1971 


Upton 


UBV photometric parallaxes 


3.25±0.20 


1972 


Golay 


Geneva photometric parallaxes 


3.30 


1972 


Iben... 


Stellar interiors 


3.23 


1973 


Kocstcr... 


Stellar interiors 


3.21 


1974 


van Altena 


Mean of secondary indicators 


3.29±0.20 


1974 


Upgren 


Trigonometric parallaxes (van Vlcck) 


3.29±0.08 


1975 


Hanson 


Compilation of methods to date 


3.19±0.15 


1975 


Klcmola... 


Trigonometric parallaxes (Lick) 


3.19±0.04 


1975 


Corbin... 


Proper motions from meridian circles 


3.42±0.20 


1975 


Hanson 


Absolute pm's wrt extragalactic 


3.18±0.16 


1977 


McAllister 


Absolute pm's corrected 


3.10±0.17 


1977 


Buchholz 


GCTSP + systematic corrections 


3.32±0.06 


1977 


Hanson 


Proper motion gradients 


3.30±0.04 


1980 


Hanson 


Weighted mean of geometric methods 


3.25±0.08 


1980 


Hanson 


Trigonometric parallaxes 


3.40±0.29 


1981 


Hauck 


Gliesc/field + Lutz-Kclkcr correction 


3.30 


1981 


Hardorp 


Masses of visual binaries 


3.47±0.05 


1982 


McClure 


Masses of visual binaries 


3.20 


1982 


Eggen 


Photoelectric photometry of 72 stars 


3.30 


1983 


Morris... 


Convergent point 


3.45±0.05 


1984 


VandcnBcrg... 


Stellar evolution theory 


3.23 


1984 


Dctwcilcr... 


Revised radial velocity 


3.26±0.11 


1985 


Cameron 


Main sequence versus Gliese stars 


3.33 


1985 


Stcfanik... 


Vrad (212 stars) -f- Hanson pm 


3.42±0.10 


1987 


Loktin... 


Proper motion geometry in FK4 


3.36±0.05 


1987 


Peterson... 


McClure data plus new photometry 


3.28±0.10 


1988 


Gunn... 


Vrad from Griffin + bulk Hanson pm 


3.35±0.07 


1988 


Heintz 


5 binaries 


3.42±0.10 


1989 


Loktin... 


Proper motion gradient 


3.37±0.07 


1990 


Schwan 


Proper motions from 44 FK5 


3.30±0.10 


1990 


Upgren... 


Parallaxes (van Vlcck, 23 stars) 


3.18±0.09 


1991 


Patterson... 


Parallaxes (McCormick, 10 stars) 


3.40±0.04 


1991 


Schwan 


Proper motions from 145 FK5/PPM 


3.45±0.06 


1992 


Morris 


Convergent point 


3.16±0.10 


1992 


Gatcwood... 


Parallax of 51 Tauri 


3.20±0.06 


1992 


Gatcwood... 


Mean parallaxes to date 


3.2 ±0.1 


1994 


Turner... 


Combined methods 


3.40±0.07 


1997a 


Torres... 


Orbital parallax 51 Tau (propagated) 


3.38±0.11 


1997b 


Torres... 


Orbital parallax 70 Tau (propagated) 


3.39±0.08 


1997c 


Torres... 


Orbital parallax 78 Tau (propagated) 


3.42±0.09 


1997a 


van Altena... 


HST FGS observations of 7 objects 


3.32±0.06 


1997b 


van Altena... 


Mean ground parallaxes to date 



Although v can in principle be determined from the ra- 
dial velocity measurements alone, its resulting direction is 
generally not well determined because of the limited angu- 
lar extent of the cluster; the usual procedure has therefore 
been to determine (v) from proper motions, and |v| from 
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radial velocities, from which A is obtained from Eq. (1) 
and 7T from Eq. (2). 

Although the method is conceptually simple, its appli- 
cation in practice is not so straightforward. Errors in the 
individual proper motions resulting from measurement er- 
rors, or defects in the proper motion system, lead to ac- 
cidental errors in ir, to an error in (v) and, ultimately, to 
a systematic bias in A depending on (b). For the Hyades, 
the streaming motion differs by only 60-70 degrees from 
that of local field stars towards the solar antapex, so that 
observational scatter in the proper motions of member 
stars, and the random motions of field stars, complicates 
membership selection based only on proper motions. 

In the basic convergent-point method it is assumed 
that the cluster is neither expanding, contracting, or ro- 
tating, that the motion of the cluster with respect to the 
field is large enough to permit accurate membership dis- 
crimination, and that the system of proper motions is in- 
crtial and without systematic errors. In his review, van Al- 
tena (1974) considered that the first two criteria were ad- 
equately satisfied, but that information on the proper mo- 
tion system was incomplete. Hanson (1975) considered the 
possibility of random motions contributing significantly to 
the stars' space velocities, as well as the effects of expan- 
sion, contraction, or rotation, concluding that any result- 
ing deviations from parallel motion are insignificant at lev- 
els affecting the distance determination by the convergent 
point method. Gunn et al. (1988) presented weak evidence 
for rotation at the levels of ^ 1 km s _1 rad -1 (projected), 
not inconsistent with these conclusions. 

The convergent point method was first applied to the 
Hyades cluster by Boss (1908), using the proper motions of 
41 suspected cluster members supplemented by three ra- 
dial velocities. The classical convergent-point method was 
further developed and discussed by Smart (1938), Brown 
(1950), and others. A systematic regression error arising 
from the quadratic form of the proper motion component 
coefficients in the normal equations, and leading to an up- 
ward revision of 7 per cent in the distance to the cluster, 
was identified by Seares (1944, 1945). 

Subsequent distance determinations using the convergent- 
point method initially appeared to be in close agreement, 
although the correspondence between the van Bucrcn 
(1952) and Wayman et al. (1965) results was later at- 
tributed in part to the use of the same proper motion 
system (van Altena 1974). Hodge & Wallerstein (1966) 
suggested that the cluster was 20 per cent farther than in- 
dicated by the proper motions - given the previous stan- 
dard distance, binary stars in the Hyades would have been 
overluminous with respect to their masses, both as com- 
pared to normal stars like the Sun, and as compared to 
models derived from stellar structure theory (Wallerstein 
& Hodge 1967). 

In the classical convergent-point method, the determi- 
nation of (v) depends only on the directions of proper 
motions, and not on their absolute values. Upton (1970) 



derived a procedure for calculating the distance directly 
from the proper motion gradients across the cluster, dis- 
pensing with the intermediate step of locating the conver- 
gent point - the cluster distance is then given by the ratio 
of the mean cluster radial velocity to the proper motion 
gradient in either coordinate. Use of this method, whose 
relevant equations can be derived by differentiating the ba- 
sic convergent-point equation (Eq. 2), has the advantage 
that more complete use is made of the proper motion data, 
while the two independently measured gradients yield two 
distance estimates whose comparison provides an indica- 
tion of the systematic and accidental errors involved. 

The accepted distance to the Hyades was revised from 
about 40 pc to about 44 pc around 1978 based on mod- 
els of the chemical composition (Koester & Weidemann 
1973), and independent astrometric and photometric re- 
sults (van Altena 1974, Hanson 1975, Eggen 1982). Hanson 
(1975) applied different formulations of the convergent- 
point method to new proper motion and cluster member- 
ship data, concluding that the errors due to different for- 
mulations of the method appeared to be quite small, with 
systematic errors in previous meridian circle proper mo- 
tions implicated as the cause of the discrepancies which 
seemed to exist between distances derived from earlier 
proper motion analyses and those resulting from a broad 
variety of other observational methods. 

Taking into account systematic magnitude effects in 
the Hanson proper motions, McAlister (1977) revised Han- 
son's distance of 48 pc downward to 43 pc, close to the 
value of 43.5 pc given by Corbin et al. (1975). Mur- 
ray & Harvey (1976) showed how all measurements of 
proper motion and radial velocity of the cluster members 
could be combined into a general solution for the cluster 
motion and the parallaxes of individual stars. A review 
of astrometric results by Hanson (1980) concluded that 
45.6 pc (distance modulus 3.30 ± 0.04) was indicated by 
the best of current data; being a weighted mean of clas- 
sical convergent-point methods and trigonometric para- 
llaxes. Meanwhile, dynamical parallaxes from visual and 
eclipsing binaries have traditionally yielded slightly higher 
--Jimcs: McClure (1982) found 49 pc, but a reanalysis of 
that and other data by Peterson & Solensky (1987) gave 
47 pc (distance modulus 3.36 ± 0.05) still slightly higher 
than the astrometric results. 

Determination of the convergent point from radial ve- 
locities was applied to the Hyades by Stefanik & Latham 
(1985), Detweiler et al. (1984), and Gunn et al. (1988), 
based on the methodology applied by Thackeray (1967) 
to Sco-Cen. 

Schwan (1990, 1991) presented the most recent deter- 
minations of the convergent point based on proper motions 
of 44 and 145 stars from the FK5 and FK5/PPM respec- 
tively, and derived a distance modulus of 3.40 ± 0.04. 

Individual trigonometric parallaxes have been pub- 
lished for certain candidate Hyades members, most re- 
cently by Patterson & Ianna (1991), by Gatewood et al. 
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Fig. 1. Distance modulus (given by m — M — — 51og7r — 5, where 7r is the parallax in arcsec) for the distance determinations 
with errors, since 1980 given in Table 1. The 'Torres' determination refers to Torres et al. (1997c). 



(1992), and in the Fourth Edition of the General Cata- 
logue of Trigonometric Stellar Parallaxes (van Altena et 
al. 1995). Recent determinations of a mean cluster dis- 
tance have been given by Turner et al. (1994), and van 
Altena et al. (1997b). High-precision orbital parallaxes 
have recently been determined for 51 Tau (vB24), 70 Tau 
(vB57), and 78 Tau (vB72) by Torres, Stefanik & Latham 
(1997a, b,c). In each case, however, their extrapolation to a 
mean cluster distance involves the use of relative ground- 
based proper motions, so that the high intrinsic accuracy 
of their orbital parallax determinations does not propa- 
gate through to a corresponding accuracy on the mean 
distance. 

A weighted mean distance modulus of 3.42 ±0.09 mag, 
based on Hubble Space Telescope FGS observations of 
seven cluster members, has been given by van Altena et 
al. (1997a). 

These recent determinations are shown, with their 
published errors, in Fig. 1. It is evident that the Hyades 
distance may still not be considered as a conventional as- 
tronomical constant. It is this uncertainty, and its atten- 
dant implications, that this paper seeks to resolve. 

The Hipparcos Catalogue provides parallaxes for all 
stars in the observing programme, of sufficient (milliarc- 
sec) accuracy not only to assign membership probabilities 
on the basis of the distances alone, but to resolve the depth 
of the cluster. The availability of annual proper motions 
with standard errors of order 1 mas yr _1 for all stars, leads 
to an opportunity to re-discuss membership on the basis 
of convergent point analysis, to probe the kinematical as- 
sumptions implicit in such analyses, and to examine, in 



combination with the parallaxes, the cluster membership 
and dynamics independently of any assumed dynamical 
model. 

3. Observational material 

3.1. Data from the Hipparcos Catalogue 

This study makes use of the final data contained in the 
Hipparcos Catalogue, which provides barycentric coordi- 
nates, incrtially-referenced proper motions, and absolute 
trigonometric parallaxes for nearly 120 000 stars (ESA 
1997). It is based on the 240 candidate Hyades mem- 
bers specifically included in the Hipparcos Input Cata- 
logue, supplemented by the astrometric and photometric 
data for all Hipparcos Catalogue objects within the range 
2 h 15 m < a < 6 h 5 m and -2° < 5 < +35° for independent 
membership studies of objects not considered as candidate 
members in the past (this region includes all previous can- 
didate cluster members). 

We stress from the outset that the Hipparcos Input 
Catalogue, on which the final Hipparcos Catalogue con- 
tents are based, is not complete to the observability limit 
of the Hipparcos observations, although specific attention 
was given during its construction to the inclusion of po- 
tentially observable candidate Hyades cluster members. 
Thus, although membership analysis can be conducted on 
previously unsuspected cluster members contained in the 
Hipparcos Catalogue, this catalogue will not contain mem- 
bers fainter than the satellite observability limit, of around 
V ~ 12 mag, nor objects omitted from the Hipparcos In- 
put Catalogue for other reasons. Clearly, an incomplete 
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Fig. 2. Hyades stars from the 'Base des Amas' (BDA). Stars contained in the Hipparcos Catalogue are displayed as filled circles 
(f90 stars). Stars not contained in the Hipparcos Catalogue, but appearing in the Tycho Catalogue, are displayed as open 
triangles (27 stars). The remaining 174 stars contained only in the BDA are displayed as open squares. 



survey will most likely result in a preferential selection of 
stars according to distance, and a biased value of the mean 
cluster distance - although it will not affect the discussion 
of the main sequence modelling when using individual dis- 
tance estimates for each object. 

The choice of Hipparcos targets in a given field of the 
sky was subject to many operational constraints which 
were, in some cases, in conflict with scientific require- 
ments. Fig. 2 illustrates the sample of Hyades stars con- 
tained in the Hipparcos Catalogue with respect to the 
global content of Hyades candidate members contained in 
the data base for stars in open clusters ('Base des Amas', 
or BDA, Mermilliod 1995), on the basis of the photomet- 
ric data, V and B — V, contained in the BDA. Stars from 
the BDA contained in the Hipparcos Catalogue arc dis- 
played as filled circles (190 stars). Stars not contained in 
the Hipparcos Catalogue, but appearing in a second major 
product of the Hipparcos mission, the Tycho Catalogue 
(ESA 1997), are displayed as open triangles (27 stars). 
The remaining 174 stars contained only in the BDA are 
displayed as open squares. The figure shows, as expected, 
the progressive incompleteness of the Hipparcos sample 
with increasing magnitude (and is also a useful demon- 
stration of the completeness of the Tycho Catalogue down 
to about V = 10.5 mag). 



Up to 114 of the 'BDA only' stars lie within 5 pc of 
the cluster centre, assuming that all lie at the mean dis- 
tance of the cluster centre. Any possible kinematical bias 
due to the present sample selection has not been studied 
in the present paper, as the distance, proper motion and 
radial velocity data for stars other than those contained in 
the Hipparcos Catalogue are either insufficiently accurate, 
incomplete, or very inhomogeneous. 

For objects contained in the Hipparcos Catalogue, 
those considered to be members or candidate members 
by one or more previous workers are listed in Table 2. 
Column a gives the Hipparcos Catalogue identifier, while 
columns b-m give the designation and membership status 
according to a number of previous workers. These are not 
the only papers where membership of particular objects 
have been discussed (see, for example, references in Grif- 
fin et al. 1988) although they represent the most substan- 
tial developments of the cluster membership studies. The 
membership status listed in columns b-m do not neces- 
sarily reflect fully the membership assignment in the orig- 
inal papers: in some cases the authors give probabilities 
for membership or include some indication of 'doubtful 
membership', which we have converted (sometimes sub- 
jectively) into a 1 or 0. Table 2 thus reflects our own un- 
derstanding of previous membership studies converted to 
a yes/no status. For full details we refer to the original 
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papers. Entries in Table 2 with '-' in columns b-m inclu- 
sive are new candidates arising from the present study, se- 
lected as described in subsequent sections (we considered 
it desirable to list all candidates sequentially in one table, 
independently of their history, and will distinguish their 
historical status by referring to these as 'previous' and 
'new' members, where appropriate). 

The previous candidates compiled in Table 2 form the 
basic list of objects for our initial studies. It is noted 
that this list contains objects already considered as non- 
members by some or even all previous workers (and which 
we will go on to confirm as non-members), while it does 
not include those objects which are considered as possible 
or secure Hyades members which are not contained within 
the Hipparcos Catalogue (as described above). Later in 
the paper, having determined the general spatial and kine- 
matical properties of the cluster on the basis of the gen- 
eral properties of the previous members, we will provide 
our own assignment of membership to this basic list (this 
final result is given in the last column of Table 2). We 
will also supplement the previous candidates by additional 
candidates selected from the Hipparcos Catalogue having 
spatial and kinematical properties in common with the 
general cluster (also included in Table 2). 

The Hipparcos Catalogue (ESA 1997) itself describes 
the details of the catalogue construction and contents, 
while recent summaries may be found in the literature 
related to the construction of the intermediate catalogue 
(Kovalevsky et al. 1995), and to the determination of the 
trigonometric parallaxes and associated errors (Perryman 
et al. 1995). 

The Hipparcos and Tycho Catalogues have been con- 
structed such that the Hipparcos reference frame coin- 
cides, to within limits set by observational uncertainties, 
with the International Celestial Reference System (ICRS), 
as recommended by the IAU Working Group on Refer- 
ence Frames (Ma et al. 1997, see also Lindegren & Ko- 
valevsky 1995). The latter system is practically defined 
by the adopted positions of several hundred extragalac- 
tic radio sources. It supersedes, although it is consistent 
with, the optical reference frame defined by the FK5 cat- 
alogue, which was formally based on the mean equator 
and dynamical equinox of J2000. The resulting deviation 
from inertial, about all three axes, is considered to be less 
than approximately 0.25 mas yr" 1 . For a discussion of the 
comparison of ground-based positions and proper motions 
with those of Hipparcos, see Lindegren et al. (1995). The 
Hipparcos Catalogue (ESA 1997, Volume 1, Section 1.5.7) 
details the relationship between the ICRS (Hipparcos) and 
J2000(FK5) frames. The epoch of the Hipparcos Cata- 
logues is J1991.25, although the provision of the full co- 
variance matrix of the astrometric solution for each star 
permits the positions, and corresponding standard errors, 
to be propagated to any epoch within the same reference 
system. 



The Hipparcos trigonometric parallaxes are absolute, 
and are considered to be free from systematic (global) er- 
rors at a level of some 0.1 mas or smaller (Arenou et al. 
1995). Studies so far suggest that the true external par- 
allax and proper motion errors are unlikely to be under- 
estimated, as compared with the formal standard errors, 
by more than about 10-20 per cent. We will demonstrate 
that the present results provide further evidence for the 
reliability of the quoted astrometric standard errors, and 
provide independent evidence for the absence of significant 
systematic errors in the parallaxes and proper motions. 
Columns (n-o) of Table 2 provide the Hipparcos parallax 
and standard errors (in mas). These values, as well as the 
remaining astrometric parameters (and correlations), are 
as published in the Hipparcos Catalogue. 

The Hipparcos Catalogue also contains detailed photo- 
metric data, including broad-band, high-precision, multi- 
epoch photometry in the calibrated Hipparcos-specific 
photometric system Hp. These are homogeneous magni- 
tudes derived exclusively from the satellite observations, 
providing the basis for detailed photometric variability 
analyses which are summarised, star-by-star, in the pub- 
lished catalogue (the Hp magnitudes were not used in the 
construction of the HR diagrams in Sect. 9 in view of the 
absence of appropriate bolometric corrections). In addi- 
tion, the catalogue includes V magnitudes, and B — V 
and V — I colour indices derived on the basis of satellite 
and/or ground-based observations. 

3.2. Binary information and radial velocity data 

Information on the binary nature of the stars in the 
Hyades is important for a variety of reasons: in addition 
to the astrophysical relevance, the confidence which can 
be placed on the kinematic or dynamical interpretation of 
the radial velocities and proper motions (and hence the 
space motions) is affected by the (known or unknown) bi- 
nary nature of the object. We have therefore attempted 
to compile the best available radial velocity and binary 
information for each object, and their inter-relationship. 

Columns p-r of Table 2 provide the radial velocity, 
standard error, and source of radial velocity, respectively. 
These columns represent the result of our literature search, 
and are supplemented by Coravel radial velocity results 
specifically acquired in the context of this study by one of 
us (JCM, column r = 24). The Coravel velocities include 
the 'standard' zero-point correction of 0.4 km s _1 (Scarfe 
et al. 1990). The Griffin et al. radial velocities given in Ta- 
ble 2 are the 'uncorrected' values given in their paper: for 
use in our subsequent kinematical studies they have been 
corrected as described in Eq. (12) of Gunn et al. (1988), 
but accounting for a sign error which is present in their 
equation (and confirmed by the authors): in their nota- 
tion, we have added a correction of —q(V) — 0.5 km s _1 
for stars fainter than V = 6 mag, and a correction of 
—0.5 km s _1 for brighter stars. In the assignment of errors 



8 



M.A.C. Perryman et al.: The Hyades: distance, structure, dynamics, and age 



Table 2. Data on the membership of the Hyades for the 282 stars in our sample, listed by various authors. Membership 
or non-membership inferred by the relevant authors are indicated by '1' or '0' in the corresponding column respectively 
(see text). Entries with '-' in columns (b-m) inclusive are new candidates proposed in this paper. Columns have 
the following meaning: (a) Hipparcos Catalogue (HIP) number; (b) van Bucren number (1952, BAN, 11, 385); 
(c) Membership according to van Bueren; (d) van Altena number (1969, AJ, 74, 2); (e) Membership according to 
van Altena; (f) Hanson number (1975, AJ, 80, 379); (g) Membership according to Hanson; (h) Pels et al. (Leiden) 
number (1975, A&A, 43, 423); van Bueren stars have the vB number + 1000; (i) Membership according to Pels et 
al.; (j) Sequential number in Table 4 of Griffin et al. (1988, AJ, 96, 172); (k) Membership according to Griffin et al.; 
(1) Schwan number (1991, A&A, 243, 386); (m) Membership according to Schwan; (n) Hipparcos parallax (mas); 
(o) Hipparcos parallax standard error (mas); (p) Radial velocity (km s _1 ); (q) Error in radial velocity (km s _1 ; # 
preceding the error indicates SB/RV (column s) with undetermined 7 velocity); (r) Source of radial velocity; (s) SB 
= spectroscopic binary, RV = radial velocity (possibly) variable; (t) H, I, M = star was previously known, or classified 
by Hipparcos, to have resolved components (from Field H56); this may overlap with the column u flag, but may also 
indicate visual or wide binary (see text for details); (u) C, G, O, V, or X = relevant part of the Hipparcos Double 
and Multiple Systems Annex, from Field H59, supplemented by S = suspected binary in Hipparcos Catalogue, from 
Field H61 (see text for details); (v) distance, d (pc), from the cluster centre defined by the 134 stars within r < 10 pc 
(see Table 3); (w) kinematic statistic c = z'S _1 z (c = 14.16 corresponding to 3cr); (x) Final membership assigned 
in this paper (0, 1); '?' indicates possible new members unclassifiable due to unknown radial velocities. 

Sources of radial velocities: (0) Radial velocity unknown; (1) Griffin et al. AJ, 96, 172 (1988); AJ, 90, 609 (1985); 
AJ, 86, 588 (1981); AJ, 83, 1114 (1978); AJ, 82, 176 (1977); A&A, 106, 221 (1982); (2) Hipparcos Input Catalogue 
(mainly from R.E. Wilson, 1953); (3) Weighted mean of rcf. 2 (39.6±1.2) and Kraft, ApJ, 142, 681 (1965, 38.4±1.5); 
(4) Kraft, ApJ, 142, 681 (1965, 37.4 ±0.4 and 36.5 ±0.5); Chcrigucnc, A&A, 13, 447 (1971, 37.3 ±0.7); (5) McClure, 
ApJ, 254, 606 (1982); (6) Torres et al., ApJ, 474, 256 (1997); (7) Mayor & Mazeh, A&A, 171, 157 (1987); (8) Kraft, 
ApJ, 142, 681 (1965); (9) Margoni et al., A&AS, 93, 545 (1992); (10) Lucy & Sweeney, AJ, 76, 544 (1971); (11) Abt 
& Levy, ApJS, 59, 229; (12) Griffin, MNRAS, 155, 1 (1971); (13) Andersen & Nordstrom, A&A, 122, 23 (1983); 
(14) Morse et al., AJ, 101, 1495 (1991); (15) Dctwciler et al., AJ, 89, 1038 (1984); (16) Weighted mean of data from 
Palmer et al., Roy. Obs. Bull., 135 (1968) and Stillwcll, PDAO, 7, 337 (1949); (17) Tomkin et al., AJ, 109, 780 (1995); 
(18) Heintz, ApJS, 46, 247 (1981); (19) Abt, ApJS, 11, 429 (1965); (20) Fekel, PASP, 92, 785 (1980); (21) Perraud, 
Journal dcs Observatcurs, 45, 361 (1962); (22) Fouts & Sandage, AJ, 91, 1189 (1986; star G83-18); (23) Strassmeier 
et al., A&AS, 72, 291 (1988); (24) New Coravel observations provided by J.C. Mcrmilliod; (25) Woolley et al., Royal 
Obs. Annals, 14, 1; (26) Hanson & Vasilevskis, AJ, 88, 844; (27) Evans, Bull. Inf. CDS, 15, 121 (1978); (28) Orbit 
recomputed by Mermilliod with period = 490 ± 1 days (from Batten) . 
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Table 2. Hyades membership compilation summary (3/4) 
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1.2 


0.14 


1 


20901 


74 


1 


504 


1 


* 




1074 


1 






14 


1 


20.33 


0.84 


+39.9 4.1 


11 


* 


* 


* 


4.1 


0.26 


1 


20916 


75 


1 


511 


1 


448 





1075 


1 


* 




104 


1 


20.58 


1.74 


+45.0 #2.5 


2 


SB 


I 




2.3 


3.93 


1 


20935 


77 


1 


536 


1 


461 


1 


1077 


1 


209 


1 


105 


1 


23.25 


1.04 


+39.90 0.11 


1 


SB 


* 





3.4 


0.77 


1 


20948 


78 


1 


544 


1 


469 


1 


1078 


1 


210 


1 


106 


1 


21.59 


1.09 


+38.62 0.24 


1 




I 


* 


1.0 


0.04 


1 


20949 


1 O 


1 




— 


* 




i n?K 
1 U i 


1 


one 


1 


57 


1 


17.08 


1.18 


+39.02 0.17 


1 


* 


* 


* 


15.2 


0.47 


1 
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7Q 

/ y 
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1 


A 7fl 
4± 1 U 




1 D7Q 

iu i y 
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107 


1 


24.19 


1.76 


+40.70 0.06 


1 


* 




* 


5.1 


0.84 


1 


20952 


* 
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474 





* 




* 




* 




7.68 


1.27 


+96.3 1.2 


2 


* 


* 


* 


83.9 


512.19 





20978 


180 





560 


1 


478 


1 


56 


1 


215 


1 


108 


1 


24.71 


1.27 


+40.97 0.06 


1 


* 


* 


* 


5.9 


1.69 


1 


20995 


80 


1 


569 


1 


481 





1080 


1 


* 




171 





22.93 


1.25 


+29.3 5.00 


18 


SB 


I 


Q 


2.9 


3.30 


1 


21008 


81 


1 




— 


* 




1081 


1 


* 




109 


1 


19.94 


0.93 


+38.0 #2.5 


2 


SB 


* 


* 


4.7 


1.28 


1 


21019 


* 




* 


— 


* 




* 




216 





* 


_ 


3.52 


1.98 


+47.89 0.26 


1 




* 


* 


237.8 


15.22 





21029 


82 


1 


584 


1 


* 




1082 


1 


* 




36 


1 


22.54 


0.77 


+41.0 1.8 


13 


* 


I 


* 


2.1 


0.75 


1 


21036 


84 


1 


591 


1 


* 


- 


1084 


1 


* 


- 


38 


1 


21.84 


0.89 


+38.8 1.2 


2 


* 


I 


* 


2.5 


0.14 


1 


21039 


83 


1 


589 


1 


493 





1083 


1 


* 


— 


37 


1 


22.55 


1.09 


+39.56 0.23 


24 


SB 


I 


* 


2.2 


0.35 


1 


21053 


OO 


1 


597 


1 
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1 
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l 
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111 


1 


24.28 


0.79 
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8 


* 
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* 
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1 
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OO 
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1 


* 


— 
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1 
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i 
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1 


22.96 


0.99 
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1 


* 


* 


* 


5.4 


1.03 


1 


21092 


























19.64 


9.61 


* * 





* 


H 




13.3 


0.95 


? 


21099 


87 


1 


* 


- 






1087 


1 


222 


1 


58 


1 


21.81 


1.25 


+40.62 0.08 


1 


* 


* 


* 


2.9 


0.38 


1 


21112 


88 


1 


625 


1 


507 


1 


1088 


1 


224 


1 


172 





19.46 


1.02 


+40.98 0.31 


1 


* 


* 


* 


5.6 


0.44 


1 


21123 






627 


1 


509 


1 


63 


1 


225 


1 


* 




23.41 


1.65 


+40.38 0.11 


1 


SB 


* 


Q 


3.8 


0.41 


1 


21137 


89 


1 


644 


1 


516 


1 


1089 


1 


* 




39 


1 


22.25 


1.14 


+36.0 #2.5 


2 


SB 


* 


* 


1.7 


2.23 


1 


21138 


191 





645 


1 


517 


1 


62 


1 


228 


1 


* 




15.11 


4.75 


+41.28 0.21 


1 


* 


* 


* 


19.9 


1.08 


1 


21152 


90 


1 


* 


— 


* 




1090 


1 


* 




143 


1 


23.13 


0.92 


+39.8 1.0 


24 


* 


* 


* 


9.4 


0.20 


1 


21179 




— 


677 


1 


532 


1 


60 


1 


* 


— 


* 




17.55 


2.97 


+41.70 #1.0 


24 


SB 


* 


* 


11.1 


2.01 


1 


21194 
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682 





541 





* 




* 




* 


_ 


9.42 


2.76 


* * 





* 


* 


* 


60.1 


69.56 





21256 


* 




* 


— 


* 




66 


1 


235 


1 


* 


- 


24.98 


1.95 


+41.39 0.20 


1 


* 


* 


* 


7.2 


1.19 


1 


21261 


* 




* 




* 




65 


1 


237 


1 
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* 


* 
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94 
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1 
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1 
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1 
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21395 






771 


1 
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* 
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1 
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1 
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* 


I 


G 
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1 


* 




* 
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1 
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SB 


I 


* 
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1 
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1 


* 
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1 
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1 


42 


1 
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0.91 


+39.86 0.29 


1 


* 


* 


* 
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1 
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1 


* 




* 
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1 
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1 


123 


1 


20.81 


1.30 


+41.86 #0.12 


1 


SB 


* 


* 


3.5 


0.74 


1 
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1 
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1 
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16 


1 
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2 


* 
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* 
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Table 2. Hyades membership compilation summary (4/4) 



HIP 


vB 




vA 




Hanson 


Pels 




Griffin 


Schwan 




# 




# 




# 




# 




# 




# 


(a) 


(b) (c) 




e) 


(f) (g) 


(h) (i) 


0) (k) 


(1) ( m ) 


21683 


108 


1 


* 


_ 


* 


_ 


1108 


1 


* 




17 1 


21684 


* 




* 


- 


* 




* 




264 





* _ 


21723 


* 




* 


— 


* 




80 


1 


266 


1 


124 1 


21741 


109 


1 




— 


* 




1109 


1 


267 


1 


60 1 


21762 


185 





" 


— 


* 




82 


1 


269 


1 


125 1 


21788 


110 


1 


* 


— 


* 




1110 


1 


270 





166 


21829 


163 





* 


— 


* 


- 


* 








* - 


21923 


* 


- 


* 


— 


* 


— 


* 


— 
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* - 


21946 


* 




* 


— 


* 




* 






u 


* - 
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* 




* 


— 


* 




O 1 


1 




u 


* - 
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111 


1 


* 




* 
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1 


* 




43 1 
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* 




* 


- 


* 




* 
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* _ 


22157 


112 


1 


* 


— 


* 




1112 


1 


* 




18 1 


22176 


164 





* 


— 


* 




* 




283 





* — 


22177 


* 




* 


— 


* 




119 


1 


285 


1 


* — 


22203 


142 


1 


* 


— 


* 




1142 


1 


284 


1 


126 1 


22221 


113 


1 




— 


* 


- 


1113 


1 


286 


1 


144 1 


22224 




— 


* 


— 


* 


— 


92 


1 




— 


127 1 


22253 


* 




* 


— 


* 




93 


1 


290 


1 


* - 


22265 




1 


* 


— 


* 




1 1 14 


1 


* 




128 1 


22271 


* 




* 


_ 


* 








291 





* _ 


22350 


115 


1 


* 


- 


* 


_ 


1115 


1 


296 


1 


61 1 


22380 


116 


1 


* 


— 


* 


_ 
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1 


298 


1 


129 1 


22394 


117 







— 


* 




* 




299 


1 


62 1 


22422 


118 


1 


* 


— 


* 




1118 


1 


300 


1 


130 1 


22446 


165 





* 


— 


* 




* 




301 





* — 


22496 


119 


1 




— 


* 


- 


1119 


1 


* 




131 1 


22505 


120 


1 




— 


* 


— 


1120 


1 


305 


1 


132 1 


22524 


121 






— 


* 




1121 








133 1 


22550 


1 


1 
1 


* 


— 


* 




1 1 n 


i 




1 
1 


134 1 


22565 


123 


1 


* 


_ 


* 




1123 


1 


* 




8 1 


22566 


143 


1 


* 


- 


* 




1143 


1 


313 


1 


135 1 


22607 


124 


1 




— 


* 




1124 


1 


* 




136 1 


22654 


* 






— 


* 


_ 


98 


1 


318 


1 


- 


22684 


145 


1 


* 


— 


* 




1145 


1 


319 





167 


22751 


125 


1 


* 


— 


* 


_ 


* 




325 





168 


22782 


146 


1 


* 


— 


* 


- 


1146 


1 


327 





174 


22805 


166 





* 


— 


* 


- 


* 


— 


328 





* — 


22850 


1 OK 


1 


* 


— 


* 




1 1 ip. 


1 


* 




137 1 


22893 


1/17 




* 


— 


* 




* 




ii i 


u 


175 


23044 


149 


1 


* 




* 


_ 


1149 


1 


336 





* _ 


23056 


148 


1 


* 


- 


* 




* 




335 





169 


23069 


127 


1 


* 


— 


* 




1127 


1 


337 


1 


138 1 


23205 
























23214 


128 


1 




— 


* 




1128 


1 


* 




44 1 


23312 


* 








* 




* 








160 


23409 


* 




* 




* 




105 


1 


341 





* _ 


23497 


129 


1 


* 




* 




1129 


1 


* 




4 1 


23498 


187 





* 




* 




107 


1 


345 


1 


139 1 


23574 


150 


1 


* 




* 




1150 


1 


346 





176 


23589 


* 




* 




* 




* 




347 





* _ 


23599 


* 




* 




* 




* 




348 





* _ 


23662 
























23701 


151 


1 


* 




* 




1151 


1 


349 
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Parallax 



Radial Velocity 



Multiplicity 



Membership 
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G 7T 


V T 


17 V 


s 


Vr H56 H59 


d (pc) 


c 


S 


(n) 


(o) 


(P) 


(q) 


« 


(s) 


(t) 


(u) 


(v) 


(w) 


(x) 


20.51 


0.82 


+35.6 


2.5 


19 


* 




* 


3.4 


2.91 


1 


9.56 


1.73 


+30.69 


0.20 


1 


* 


* 


* 


58.7 


63.17 





23.95 


1.63 


+42.50 


0.19 


1 


* 


* 


* 


5.9 


0.53 


1 


15.96 


1.36 


+41.34 


0.16 


1 


* 


* 


* 


17.7 


0.52 


1 


23.65 


2.53 


+40.90 


0.17 


1 


SB 


M 


C 


4.7 


0.76 


1 


19.48 


1.26 


+35.85 


0.05 


1 








8.3 


7.85 


1 


5.88 


0.97 


+35.7 


1.2 


2 


* 


* 


* 


126.7 


47.56 





23.23 


1.25 


+ 13.57 


0.20 


1 


* 




* 


4.6 


211.77 





21.10 


2.22 


+ 14.49 


0.23 


1 


SB 


* 


* 


3.5 


178.36 





21.48 


1.84 


+24.47 


0.34 


1 


* 


* 


* 


5.4 


82.06 





20.73 


0.88 


+39.6 


0.5 


2 


* 




* 


5.9 


1.18 


1 


9.08 


1.79 


+25.94 


0.11 


1 


* 


* 


* 


64.1 


71.72 





12.24 


0.86 


+43.0 


1.0 


11 


SB 




* 


36.1 


4.32 


1 


10.81 


0.94 


+44.11 


0.10 


1 


SB 


* 


* 


46.5 


62.75 





22.45 


2.32 


+43.16 


0.25 


1 


* 


* 


* 


11.1 


0.36 


1 


19.42 


1.09 


+42.42 #0.71 


1 


SB 


* 


* 


6.5 


0.76 


1 


26.26 


1.04 


+42.47 #0.11 


1 


SB 


* 


G 


10.5 


5.44 


1 


24.11 


1.72 


+40.32 #0.09 


1 


SB 


* 


* 


6.0 


0.38 


1 


15.74 


1.98 


+41.78 


0.23 


1 


* 


* 


* 


18.7 


1.09 


1 


19.81 


1.43 


+39.8 


#0.4 


15 


SB 


* 


* 


5.9 


1.79 


1 


22.07 


2.03 


+40.30 


0.17 


1 


* 


* 


* 


8.5 


4.18 


1 


19.30 


1.67 


+41.84 #0.44 


1 


SB 


* 


G 


7.9 


0.72 


1 


21.38 


1.46 


+41.62 


0.15 


1 


* 


* 


* 


4.4 


0.97 


1 


18.96 


1.62 


+40.60 


0.31 


1 


SB 


* 


* 


10.4 


0.34 


1 


19.68 


0.96 


+42.04 


0.14 


1 


* 


* 


* 


6.3 


0.53 


1 


13.26 


1.11 


+31.84 


0.14 


1 


* 


* 


* 


30.3 


35.59 





22.96 


1.17 


+41.40 


0.16 


24 


SB 


* 


G 


5.1 


1.31 


1 


23.64 


0.99 


+42.34 #0.33 


1 


SB 


* 


g 


6.0 


1.75 


1 


19.30 


0.95 


+42.74 


0.17 


1 


SB 


* 


* 


7.2 


0.80 


1 


20.15 


1.14 


+42.44 #0.17 


1 


SB 


I 


C 


7.4 


0.83 


1 


17.27 


0.82 


+36.8 


1.2 


2 


* 


I 


* 


12.8 


4.10 


1 


17.14 


1.00 


+42.92 


0.19 


1 


* 


* 


* 


13.1 


1.19 


1 


23.91 


1.04 


+39.83 


0.24 


1 


SB 


I 




6.7 


1.87 


1 


18.93 


2.02 


+42.88 


0.25 


1 


* 


* 


* 


8.4 


0.58 


1 


12.14 


2.22 


+48.53 


0.08 


1 


* 


I 


c 


36.7 


17.34 





11.62 


1.95 


+48.72 


0.18 


1 


* 


* 


* 


40.9 


16.84 





14.82 


0.88 


+57.10 


0.23 


1 


* 


* 


* 


22.0 


95.23 





5.52 


1.29 


+ 19.14 


0.33 


1 


SB 


* 


* 


135.2 


114.49 





14.67 


0.95 


+38.4 


2.0 


8 


* 


* 


* 


22.9 


1.27 


1 


9.66 


1.43 


-30.57 


0.26 


1 


* 


* 


* 


58.3 


993.75 





12.62 


1.89 


+37.96 


0.17 


1 


* 






37.6 


13.12 


1 


14.29 


1.48 


+60.93 


0.17 


1 


* 


* 


* 


26.8 


109.89 





19.66 


1.62 


+43.68 


0.16 


1 


* 


* 


* 


7.9 


0.80 


1 


10.73 


1.66 


* 







* 


I 


c 


50.5 


8.53 


7 


23.09 


0.83 


+42.5 


1.5 


24 


* 


I 


* 


6.7 


0.30 


i 


16.77 


1.79 


+42.21 


0.40 


24 


* 


* 


* 


18.7 


1.18 


1 


11.39 


1.66 


+ 78.64 


0.36 


1 


* 


* 


* 


45.5 


358.75 





20.01 


0.91 


+38.0 


1.7 


8 


* 


* 


* 


8.9 


1.29 


1 


18.44 


1.66 


+43.51 


0.19 


1 


* 


* 


G 


11.1 


0.08 


1 


2.26 


1.27 


+30.22 


0.19 


1 


* 


* 


* 


396.8 


31.53 





5.30 


0.81 


+49.35 


0.21 


1 


* 


* 


* 


143.2 


40.64 





3.98 


1.61 


+ 111.51 
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from the compilation of Griffin et al. (1988) we adopted 
the 'internal error' quoted in their Table IV for stars with 
3 measurements or less. 

Column s indicates whether the object has been clas- 
sified as a spectroscopic binary (SB), or as (possibly) vari- 
able in measured radial velocity and therefore indicative 
of a possible spectroscopic binary (RV), according to the 
given source. Radial velocities are systemic (7) velocities 
where available. If the radial velocity has been noted as 
variable, or if no 7 velocity is available, a precedes the 
radial velocity error, indicating that it should be viewed 
with caution for dynamical studies of the cluster. Of the 
previous candidates in Table 2, 71 are classified as spectro- 
scopic binaries, 37 of which have a 7 velocity determined; 
these are mostly from the work of Griffin et al. (1988, and 
references therein), some are from refs. 10-11 accompany- 
ing Table 2, with others as referenced individually in the 
key to Table 2. [To assist cross-referencing to the results of 
Griffin et al., column j of Table 2 uses the sequential num- 
ber of the object in Table IV of Griffin et al.; this sequen- 
tial numbering takes account of two 'blocks' containing 6 
objects (rather than the usual 5) in their Table IV]. 

Columns t-u provide information on the (possible) bi- 
nary nature of the star taken from the Hipparcos Cat- 
alogue. Each entry in the Hipparcos Catalogue includes 
duplicity /multiplicity information derived from the obser- 
vations: very broadly, systems are resolved if their sepa- 
rations are larger than approximately 0.1 arcsec and their 
magnitude differences smaller than about 3 mag. For such 
systems the catalogue provides detailed information on the 
components in Part C of the Double and Multiple Sys- 
tems Annex. Additional information, such as suspected 
duplicity inferred from the astrometric residuals, or ob- 
served photocentric acceleration implying the presence of 
short-period orbital systems, has been derived from the 
observations, and relevant catalogue entries are flagged 
accordingly and assigned to distinct parts of the Double 
and Multiple Systems Annex. 

Column t is taken from Field H56 of the Hipparcos 
Catalogue and indicates that a CCDM identifier, denot- 
ing entries in the 'Catalogue of Components of Double 
and Multiple Stars' (Dommanget & Nys 1994), has been 
assigned to the catalogue entry: H indicates that the sys- 
tem was determined as double or multiple by Hipparcos 
(previously unknown); I that the system was identified as 
double in the Hipparcos Input Catalogue; and M that the 
system had been previously identified as double, but not 
recorded as such in the Hipparcos Input Catalogue. Col- 
umn u is taken from Field H59 of Hipparcos Catalogue, 
and indicates which part of the Double and Multiple Sys- 
tems Annex the entry has been assigned to: 'C indicates 
components are resolved, 'G' indicates that a non-linear 
motion of the photocentre has been detected, 'O' that the 
entry is classified as an orbital system, 'V that the entry 
is inferred to be double from a correlation between photo- 
centric motion and photometric variability, and 'X' that 



the entry is likely to be an unclassified (close) double or 
multiple system. 'S' in this column is taken from Field H61 
of the Hipparcos Catalogue and indicates, somewhat inde- 
pendently, that the entry is a suspected binary. Informa- 
tion in column t may overlap with that given in column u, 
i.e. it may simply indicate that detailed information on 
components is given in Part C of the Double and Multiple 
Systems Annex, but it may also indicate that the star is 
a component of a visual or wide binary (and not included 
in Part C). 

In summary, 'C in column u indicates that compo- 
nents of a double or multiple system have been resolved 
by the Hipparcos observations, while 'G', 'O', 'V, 'X', or 
'S' indicates that the star is, or may be, a close binary sys- 
tem. The implications for the measured radial velocities 
will be taken into account in the discussions of member- 
ship and dynamics of the cluster. 



4. Proper motions and the convergent point 

Before developing the analysis of the 6-dimensional (po- 
sition and velocity) data set provided by the Hipparcos 
proper motions and parallaxes in combination with the 
ground-based radial velocities, it is instructive to refer to 
the most recent determinations of the convergent point 
based upon the best-available ground-based data, and to 
examine the sensitivity of the resulting analysis to the ac- 
curacy of the available proper motion data. We have not 
investigated all numerical implementations of the conver- 
gent point method using the Hipparcos data - the ob- 
jective in this section is merely to gain insight into the 
performance and consistency of the classical convergent 
point methods. 

Figs 3(a) and (b) assemble the set of stars selected by 
Schwan (1991) for his convergent point analysis, and show 
the motions of the selected stars - including the region of 
the resulting convergent point - based on Schwan's data 
(a), and for the Hipparcos Catalogue data for the same 
selection of stars (b). Note that Schwan's data are referred 
to B1950(FK5). 

Inspection of the less well-defined convergent point 
apparent in Fig. 3(b) compared with that of Fig. 3(a) 
could lead to the erroneous conclusion that the Hipparcos 
proper motions are of a degraded accuracy compared with 
those used by Schwan (1991). The correct explanation is, 
rather, that for any given set of proper motions and as- 
sociated errors, the convergent point analysis selects as 
candidate cluster members those having a minimum dis- 
persion about the selected convergent point. Evidently, 
for a given set of candidate members which have been se- 
lected according to a given, but erroneous, distribution 
of proper motions, an improvement in the correspond- 
ing proper motion accuracies will not necessarily result 
in a 'tightening' of the previously-determined convergent 
point for the same selection of candidate members. That 
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a 0, FK5, B1950) a (°, ICRS, J1991.25) a (°, ICRS, J1991.25) 

Fig. 3. (a) The set of stars selected by Schwan (1991) for his convergent point analysis, and according to his adopted membership 
criteria, showing the positions (solid circles) and motions of the selected stars on their great circles (lines) - including the region 
of the resulting convergent point (intersection of the lines) - based on Schwan's data; (b) shows the Hipparcos Catalogue data 
for precisely the same selection of stars (note the difference in reference systems and epochs); (c) shows the stars selected from 
the Hipparcos Catalogue data according to Jones' (1971) method. Notice the much 'tighter' distribution of stellar motions 
compared with (b). 



the Hipparcos Catalogue data result in an increase in the 
scatter for the same selection of stars is a direct conse- 
quence of retaining a sub-optimum sample of stars on the 
basis of their (imprecise) proper motions. A revised analy- 
sis, as we will demonstrate, leads to a different convergent 
point, and a correspondingly different selection of stars. 
The consequences for the determination of the individual 
parallaxes of the candidate cluster members, and the re- 
sulting mean cluster distance, then follow directly from 
Eqs. (1) and (2). 



Our implementation of the convergent point method 
applied to the Hipparcos proper motions used the maximum- 
likelihood technique described by Jones (1971), which has 
the merit of locating the convergent point and, simulta- 
neously, the corresponding cluster members. In essence, 
members are searched for amongst the set of stars show- 
ing the clearest converging motions. Applying it to the 
Hipparcos proper motions of the previous candidate mem- 
bers listed in Table 2 (those with entries in columns b-m 
which, we recall, is our basic starting point containing a 
large proportion of possible cluster members) resulted in 
about one half of the stars (113) being selected as cluster 
members in the first iteration. Two further applications 
of the same method to those stars not selected in the first 
iteration provided 17 and 18 additional candidate mem- 
bers in the second and third steps respectively. As shown 
in Fig. 3(c) the convergence of stars selected in the first 
step is now much tighter than seen in Fig. 3(b), although 
the successively selected groups have different convergent 
points: at (a, S) = (98?6, 6?4), (95?1, 8?3), (96?6, 5?8), re- 
spectively (ICRS, epoch J1991.25). 
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Fig. 4. The differences in the proper motions components 
(HIP - FK5 etc.) for the membership candidates selected by 
Schwan (1991), normalised to the combined standard errors in 
each component. Proper motion components given by Schwan 
(1991) in B1950(FK5) coordinates have been transformed to 
J2000(FK5) for comparison with the ICRS (Hipparcos) data. 
The good agreement between the FK5 and Hipparcos proper 
motions, with degraded accuracies from Schwan's candidates 
with proper motions only from N30 or PPM, is evident. 



The explanation for the absence of a unique convergent 
point is that the methods of Jones and Schwan for select- 
ing members of the cluster (and convergent point methods 
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in general) rely on finding stars which show the least de- 
viations of the position angle of their proper motion, 9, 
with respect to the direction from the stellar position to 
the convergent point, 9 C . In reality each star has its own 
'convergent point', given by the direction of its space mo- 
tion, and the difference A9 = 9 — 9 C will be small if the 
stellar convergent point is close to the great circle con- 
necting the cluster centre to the cluster convergent point. 
As a result, the convergent point method tends to select 
stars with space motions that lie in the plane defined by 
the great circle passing through the cluster centre and the 
convergent point. 



H -'K"i/KK4SUP: 

□ K30 

□ PPM 



^(/O (Schwan, mas/yr) 



(Schwan, mas/yr) 



2 4 

ct(^J (HIP. mas/yr) 



<j(m 4 ) (HIP, mas/yr) 



Fig. 5. The distribution of standard errors in the proper mo- 
tion components for the membership candidates selected by 
Schwan (1991). Top: as used by Schwan (1991), derived from 
FK5, N30, and PPM. Bottom: the distribution of standard 
errors for the same objects from the Hipparcos Catalogue, ac- 
cording to their appearance in the previous source catalogues. 
The slightly better Hipparcos proper motion accuracies for the 
subsets of the FK5 and N30 objects reflects their generally 
brighter apparent magnitudes. 

If there is a significant velocity dispersion in the clus- 
ter the convergent point membership selection will lead 
to an artificial flattening of the distribution of candidate 
members in velocity space, which may in turn lead, for 
example, to spurious inferences of rotation. Conversely if 
there is significant systematic structure in the internal ve- 
locities of a cluster, the convergent point method may lead 
to a spatial bias in the selection of candidates. This would 
happen if the cluster possessed a significant component 
of rotation with the extreme internal velocities located 
primarily in a plane perpendicular to the great circle con- 
necting the cluster centre and convergent point. 

Schwan's proper motions were from mixed sources, 
drawn from the FK5/FK4Sup, N30, and PPM catalogues. 



That the differences in membership selection resulting 
from Schwan's values and the present Hipparcos values are 
arising from the different quality of the available proper 
motions is evident from Fig. 4, which illustrates the differ- 
ences in the proper motion components (HIP— FK5 etc) 
for the membership candidates selected by Schwan (1991), 
normalised to the combined standard errors in each com- 
ponent. The corresponding distributions of the standard 
errors are shown in Fig. 5. The generally very good agree- 
ment between the FK5 and Hipparcos proper motions, 
with degraded accuracies for Schwan's candidates with 
proper motions only from N30 or PPM, is evident. It 
should be borne in mind that the differences between, 
for example, FK5 and Hipparcos proper motions may 
partially reflect true differences in the measured proper 
motions of astrometric binaries where the FK5 proper 
motions reflect the long-term photocentric motion, with 
the Hipparcos measurements made over a period of only 
3.5 years carrying information on orbital perturbations 
over these time scales. Indeed, there is some evidence that 
the proper motion differences between Hipparcos and FK5 
exceed their combined standard errors, at least in a sta- 
tistical sense (Wielen 1997). 

We have identified that heterogeneous ground-based 
proper motions will affect the determination of the conver- 
gent point, the membership determination, and hence the 
cluster distance modulus. Although, alone, a systematic 
error of approximately 1 mas yr _1 in the proper motion 
system would be needed (Eq. 2) to account for distance 
errors of 1 per cent, the sin A term in Eq. (2) results in 
a greater sensitivity of the distance estimate to a com- 
bination of the proper motions and adopted convergent 
point. We will return to a discussion of the convergent 
point based on the Hipparcos data, and the consistency 
between the distances inferred from the convergent point 
(derived from the ground-based and Hipparcos proper mo- 
tions) and the Hipparcos-based trigonometric parallaxes, 
in Sect. 6.1. 



5. Membership determination and mean cluster 
distance 

5.1. Determination of positions and space motions 

The discussions of Sect. 4 also illustrate the point acutely 
evident to previous workers that membership selection 
based on proper motion data alone, however accurate, may 
also lead to erroneous inferences about cluster member- 
ship if significant departures from strictly parallel motion 
exist within the cluster. 

With the availability of the full 6-dimensional position- 
velocity data, based on the 5 astrometric parameters pro- 
vided by Hipparcos supplemented by the stellar radial ve- 
locity when available, we are in a position to examine 
membership based on stricter spatial and kinematic cri- 
teria than has been possible hitherto. 
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We first assemble the equations used for the transfor- 
mation between equatorial and Galactic coordinates, and 
for the determination of space velocities based on the ob- 
served proper motions, parallaxes and radial velocities. 
The transformation between the equatorial and Galactic 
systems is given by: 



[x G yc zg] = [x y z] A g 



(3) 



where [x y z] and [xq yc z g] are the basis vectors in the 
equatorial and Galactic systems respectively (x is the unit 
vector towards (a, 6) = (0, 0), y is the unit vector towards 
(+90°, 0), and z the unit vector towards S = +90°), and 
where the matrix Ag relates to the definition of the Galac- 
tic pole and centre in the ICRS system. Currently no def- 
inition of this relation has been sanctioned by the IAU, 
and we adopt the following definition proposed by the 
Hipparcos project (ESA 1997, Volume 1, Section 1.5.3), 
using as celestial coordinates of the north Galactic pole in 
the ICRS system: 



a G = 192?859 48 
5 G = +27?128 25 



(4) 



with the origin of Galactic longitude defined by the Galac- 
tic longitude of the ascending node of the Galactic plane 
on the equator of ICRS, taken to be: 



In = 32?93192 



(5) 



Eqs (4) and (5) preserve consistency with the previous 
B1950 definition of Galactic coordinates (Blaauw et al. 
1960) to a level set by the quality of optical reference 
frames prior to Hipparcos, accounting for both the trans- 
formation to the J2000(FK5) system (cf. Eq. (33) of 
Murray 1989) and then to the ICRS (Hipparcos) system 
by application of the orientation difference between the 
Hipparcos and FK5 Catalogues. 

These values of the angles a Gl S G and Iq are to be 
regarded as exact quantities. From them, the transforma- 
tion matrix Aq may be computed to any desired accuracy. 
To 8 decimals the result is: 



A G 



-0.054 875 56 
-0.873 43709 
-0.483 835 02 



+0.494109 43 
-0.444 829 63 
+0.746 982 24 



the astronomical unit expressed in the appropriate form; 
7T is the parallax expressed in mas, and 

/^ck* — 1~^ol cos 

and [is are the proper motion components expressed in 
mas yr _1 . The Doppler factor, k = (1 — Vr/c) _1 , is re- 
quired to account rigorously for light-time effects in the 
calculation of the space velocity in terms of the observed 
proper motion and radial velocity. 

In the equatorial system the components of the normal 
triad [p, q, r] are given by the matrix: 



R 



-0.86766615 
-0.198 076 37 ] (6) 
+0.455 983 78 



If b denotes the barycentric position of the star, measured 
in parsec, and v its barycentric space velocity, measured 
in km s _1 , then: 



b = A p u/ir 
and: 

v = (p(i a *A v /ir + q[i$A v /Tr + rV R )k 



(7) 



(8) 



where u is the unit vector in the barycentric direction, 
[p q r] is the normal triad defined below, and A p = 
1000 mas pc and A v = 4.74047... km yr s" 1 designate 




sin S cos a cos S cos a 
sin 5 sin a cos S sin a 
cos S sin S 



(9) 



The equatorial components of b and v may thus be writ- 



ten: 




(10) 



and: 



kfl a ^A v 1 7T 

R | kfisA v /-K 
kV R 



(11) 



The Galactic components of b and v are obtained through 
pre-multiplication by A G . In the following we ignore the 
Doppler correction factor, and set k = 1. 

5.2. Preliminary membership selection 

At this point, our notion of cluster 'membership' is inten- 
tionally vague, based only on some general preconceptions 
about the uniformity of the space velocities in the central 
region. As we shall see, realistic A^-body simulations pre- 
dict, or reflect, dynamical properties such as mass segre- 
gation and cluster evaporation as members diffuse beyond 
the cluster tidal radius, or are ejected in dynamical inter- 
actions closer to the cluster core. Thus we might expect an 
[increasing dispersion of the space velocities with increas- 
ing distance from the cluster centre. 

Our approach will therefore involve the following steps: 
(i) assign preliminary membership based on rather non- 
rigorous spatial and kinematical criteria; (ii) estimate a 
preliminary centre of mass and centre of mass motion; 
(iii) examine the displacements and velocity residuals of 
each candidate member with respect to these preliminary 
reference values; and finally (iv) refine the membership 
criteria accordingly, once the preliminary spatial and ve- 
locity structure becomes more evident. Unlike previous 
implementations of the moving cluster method, we need 
not assume anything about the degree to which the cluster 
members participate in uniform parallel motion in space; 
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rather, we will be able to examine directly the assump- 
tions on which these methods have been invoked and, in 
particular, whether there is evidence for cluster rotation, 
expansion, or shear. 

Our starting list of possible candidate members is 
given in Table 2, which includes all objects in the Hipparcos| 
Catalogue which have, in the past, been considered as 
(possible) Hyades members or assigned a reference num- 
ber in the quoted sources. Preliminary membership was 
assigned to a subset of these previous candidate members 
based on approximate limits placed on the parallax, the 
radial velocity, the object's position in the proper motion 
vector point diagram, and distributions of the Galactic 
components of b and v. This resulted in the elimination 
of the most obvious non-members and a list of 188 prelim- 
inary members. The subsequent stages of the membership 
selection are insensitive to this preliminary membership, 
so neither the details of this selection process, nor the list 
of preliminary members, are provided. 

5.3. Preliminary determination of the centre of mass 

The preliminary members show a projected spatial dis- 
tribution extending over 10-20 pc in each coordinate on 
the sky; thus for a distance of ~ 45 pc, and for median 
Hipparcos parallax errors of o~ n ~ 1 mas (~ 1.5 pc at this 
distance), the depth of the cluster is clearly resolved by 
the Hipparcos parallaxes. Our next step is to determine 
an approximate centre of mass for the cluster, not with 
the ultimate goal of determining a mean cluster distance, 
but rather in order to establish a well-defined cluster refer- 
ence point and mean space motion to which more careful 
membership assignment may be referred. 

In order to determine the centre of mass of the clus- 
ter, masses of single stars have been determined using 
a reference isochrone from the models of Schaller et al. 

(1992) . The point on the isochrone nearest to the observed 
values of My and B — V was determined for each star, 
and a mass assigned corresponding to that point on the 
isochrone. The masses of spectroscopic binaries were taken 
from the literature or, if no published mass was avail- 
able, an estimate of the minimum mass was made based 
on the luminosity of the primary. For resolved binaries, 
masses of the secondary were assigned according to the 
mass-luminosity calibration given by Henry & McCarthy 

(1993) . For components of binaries with separations larger 
than about 20 arcsec, the number density criterion given 
in Brosche et al. (1992) was used to infer whether they 
are physically associated with the primary or not. Opti- 
cal companions cross-referenced in the Hipparcos Input 
Catalogue through their CCDM identifier (see Sect. 3.2), 
mostly at separations of more than 100 arcsec, were not 
considered to be members of the Hyades, and were not 
included in the determination of the mass of the system 
- although the proper motions of these components may 
be listed in the Hipparcos Input Catalogue as identical 



to those of the 'primary', it is more probable that most 
of these are background objects. For an estimation of the 
errors on the position of the centre of mass, and for the 
dynamical investigations in Sect. 8, associated standard 
errors were arbitrarily assigned to be 0.1 M Q for single 
stars and 0.5 M for double stars. 

The centre of mass was then determined as Smibi/Emi.| 
To avoid outliers in the space positions affecting the centre 
of mass, stars located in the central regions of the cluster, 
defined (in parsecs and Galactic coordinates) by: 

-50 <b x < -30 

-10 <b y < +10 (12) 

-25 <b z < -10 

were selected for the determination of the centre of mass. 
Of the 188 preliminary members 142 lie in this central re- 
gion. The resulting (preliminary) centre of mass in Galac- 
tic coordinates (in pc) is shown in the first line of Table 3. 

The same 142 stars (of which 141 have a measured ra- 
dial velocity) were used to derive the centre of mass mo- 
tion. The derived velocity components in Galactic coordi- 
nates, and total space motion, are also given in the first 
line of Table 3. Assigning the binaries half the weight of 
single stars, to account for the larger uncertainties in their 
space motions, or using the inverse of the standard errors 
as weights, results in coordinates of the centre of mass 
which differ by no more than 0.3 pc in each component 
from the unweighted results, and in a mean velocity within 
a few tenths of km s _1 . We conclude that these results are 
rather insensitive to the weighting scheme adopted. 

5.4- Final membership selection 

The space velocity derived in the previous subsection can 
now be used to refine the membership criteria, and to de- 
termine additional candidate Hyades members (not pre- 
viously considered as members according to columns b-m 
of Table 2) from the field around the cluster according to 
kinematic criteria. As stated in Sect. 2 the expected in- 
trinsic velocity dispersion around the mean cluster motion 
is less than 1 km s _1 ; a dispersion of around 0.2 km s _1 
would be expected for a Plummer potential with a core ra- 
dius of 4° (Gunn et al. 1988 derived a corresponding core 
radius of 3.15 pc) and a mass of about 400 M Q . For more 
realistic, simple stellar systems, a larger dispersion, of or- 
der 0.4 km s _1 for a half-mass radius of 5 pc (see Sect. 8) 
may apply (Binney & Tremaine 1987, Eq. 4-80b). In se- 
lecting cluster members, at least this intrinsic dispersion 
must be considered. But numerous other effects may be 
responsible for a dispersion of the velocities around the 
mean cluster motion. In addition to observational errors, 
including the contribution of different zero points for dif- 
ferent radial velocity sources, the presence of (undetected) 
binaries will inflate the intrinsic velocity dispersion due to 
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their orbital motion, affecting both the observed proper 
motion and radial velocity distribution. 

The combined effects are clearly observed in the ve- 
locity distributions of the preliminary members. For the 
142 candidate members located in the central region of the 
cluster (of which 141 have a measured radial velocity com- 
piled in Table 2), the standard deviation in the Galactic 
components of v are 2.3, 1.9 and 2.1 km s _1 . The me- 
dian errors in the three components are 0.59, 1.04, and 
0.84 km s _1 , suggesting that the internal dispersion is re- 
solved. However, taking only the 51 single stars that have 
a radial velocity measured by Griffin et al. 1988 (essen- 
tially providing a single source of radial velocities with a 
carefully defined origin) the standard deviations are found 
to be 0.4, 0.8, and 1.0 km s _1 (in this case we have cho- 
sen half the inter-quartile range, which provides a more 
robust estimator for small samples of data), with median 
errors of 0.48, 1.29, and 0.92 km s _1 . Thus the internal 
velocity dispersion is not resolved significantly in any of 
the components of v. 

Hence we assume at this stage that all cluster members 
move with the same velocity vector even though, due to 
the use of different radial velocity sources and the presence 
of binaries, the velocity data do show considerable disper- 
sion. This has been accounted for by assigning an uncer- 
tainty to the relative velocity of each star with respect to 
the mean cluster space motion, which is the sum of the 
standard errors and covariances in the centre of mass mo- 
tion and the observed standard errors and covariances in 
the individual data. We stress that this uncertainty does 
not reflect the physical internal dispersion of the velocities 
in the Hyades, but rather the quality of the data presently 
available. 

The membership selection then proceeds as follows. 
For each star we can calculate, from Eq. (11), the expected 
values of the transverse and radial velocities for a star 
moving with the common cluster motion: 



and: 



V a * 

Vr 



R' 



(13) 



o \ " z / c 

where V a * = n a *A v /it, V$ — [isA v /ir, and Vr are the ve- 
locity components in equatorial coordinates, and the sub- 
script C refers to the centre of mass motion. Similarly, 
from the observed values of ir, fjL ait , n$, and Vr we can cal- 
culate the values of the observed transverse and radial 
motions (V^*, Vg, Vr). 

A comparison between the expected and observed ve- 
locities requires an evaluation of the associated covariance 
matrices. If the covariance matrix associated with vector 
x is C x , then the covariance matrix of y = F(x) is given 
by C y = JC X J' where J is the Jacobian matrix associated 
with the transformation from x to y. The Jacobian matri- 
ces for the transformations from vc to (V a *, V$, Vr)o, and 
from {n,fj, a *,iis,V R ) to {V a *,Vs,V R ) are: 

J = R' (14) 



/vr 2 A v /ir 
-fi s A v /ir 2 





(15) 



respectively. We will then consider a star to be a candidate 
Hyades member based on these kinematic criteria if the 
difference between the expected and observed velocities 
lies within a certain combined confidence region of the 
two calculated vectors. Assuming the two are statistically 
independent, the combined confidence region is described 
by the sum of the two covariance matrices, S, as: 



c = z'XT 



(16) 



where z is the difference vector, and c is dimcnsionlcss. 
In order to account more easily for objects for which 
the radial velocity was unavailable, the membership se- 
lection was carried out in equatorial coordinates, with 
vc = (—6.22,44.97,5.36) km s _1 , this value correspond- 
ing to the velocity in Galactic coordinates given in the first 
line of Table 3. The uncertainty in the centre of mass mo- 
tion corresponding to the 142 candidate members derived 
in Sect. 5.2 is: 



+2.40 -0.18 +0.04' 
-0.18 +2.45 +0.17 
+0.04 +0.17 +1.26 



(17) 



where the diagonal elements are the standard errors in 
km s _1 and the off-diagonal elements are the associated 
correlation coefficients. For stars with unknown radial ve- 
locity the same procedure, excluding the component Vr, 
can be applied. In this case, we use information restricted 
to the tangential velocity components, using only the 2x3 
sub-matrix of Eq. (15). The quantity c is distributed ac- 
cording to a \ 2 distribution with 2 or 3 degrees of freedom, 
depending on the dimensions of z. An adopted 99.73 per 
cent confidence region (corresponding, somewhat arbitrar- 
ily, to ±3cr for a one-dimensional Gaussian) corresponds 
to a value of c = 14.16 for the full 3-dimensional differ- 
ence vector (P — 0.9973 for x 2 = 14.16 and v = 3), or 
to c = 11.83 for the 2-dimensional difference vector if the 
radial velocity is unknown (P = 0.9973 for \ 2 = 11.83 
and v = 2). We note that the uncertainties on the mean 
cluster motion are correlated, and return to this point in 
our discussion of the velocity distribution of the cluster 
members. 

Applying the selection procedure to the 5499 stars in 
the Hipparcos Catalogue located in the area of sky noted 
in Sect. 3, and using Vr from Table 2 when available, or 
from the Hipparcos Input Catalogue compilation, when 
available, results in a list of 1027 candidate members. 
Fig. 6 shows the distribution of tt for these stars - the 
Hyades cluster corresponds to the peak in the parallax 
distribution between 15 and 27 mas. Most of the 'new 
candidates' from this list of 1027 stars have small pa- 
rallaxes with relatively large ov/V, and small expected 
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Fig. 6. The distribution of tt for the 1027 stars described in the 
text. A subsequent membership selection criterion of ir > 10 
mas was applied to these candidates. 

proper motions, due to a combination of their large dis- 
tances and/or their location close to the cluster's conver- 
gent point. Radial velocity information is largely absent 
for these additional objects, which are generally likely to 
be unassociated with the Hyades cluster, and identified 
here as 'possible' candidates (at least we cannot exclude 
them on the basis of our data) simply due to the large 
uncertainties in the observational material. 

The distribution of n offers no unambiguous criteria for 
further constraining possible membership. In the follow- 
ing we simply exclude objects with tt < 10 mas from fur- 
ther consideration, leaving a list of 218 candidate Hyades 
members; adopting a membership threshold of tt > 8 mas 
would result in 246 candidate members. Of the 218, 179 
are in the list of previous candidate members from Table 2, 
with the 39 'new' candidates also listed in Table 2 ('-' in 
all of columns b-m). The remaining 64 out of the 282 
objects in Table 2 are considered now as non-members. 

In compiling the list of new candidates in Table 2 we 
have already made use of additional radial velocities ac- 
quired by one of us ( JCM) using Coravel, which were valu- 
able in excluding some of the stars located at large dis- 
tances as members - while Griffin et al. (1988) estimated a 
1 per cent probability that a field star will exhibit a radial 
velocity indistinguishable from that of a cluster member, 
the combination of radial velocity and proper motion data 
provides an almost unambiguous membership criterion for 
objects participating in the same overall space motion. 
These new radial velocities have already been introduced 
into Table 2, as indicated in the notes to the table. New 
velocities acquired during late 1996 were already used to 
suppress intermediate candidates: this is the case for HIP 
10920, 11815, 15288, 15406, 19757, 22802, 22809, 23810, 
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Fig. 7. The HR diagram 


represented by the 1027 stars selected 



according to the kincmatical criteria (all points), and those 
retained according to the additional criterion ir > 10 mas (•). 



25419, while a further 23 stars with tt < 10 mas could 
also be excluded as candidate members based on the new 
radial velocities. The final status is that all of the previous 
179 candidate members in Table 2 now have a known ra- 
dial velocity, while 18 of the 39 'new' candidates in Table 2 
have a known radial velocity. More definitive membership 
assignment would clearly benefit from the acquisition of 
radial velocities for those candidates for which the full 3-d 
space motion is presently unknown. 

The process of recalculating the centre of mass af- 
ter each membership selection step, and further refin- 
ing it once new radial velocities were acquired and pu- 
tative members rejected, leads to velocity changes by only 
~ 0.1 km s _1 or so in each component. The resulting sen- 
sitivity of membership selection to changes in the result- 
ing systemic space motion is also small: varying all three 
velocity components by ±0.5 km s _1 (over a 3 x 3 x 3 
grid), the number of selected stars in the inner 10 pc only 
changes by 3 or 4 at most (out of 134). Taking steps of 
1 km s _1 the number of stars changes by ±5. So we may 
be confident that the precise value of the adopted space 
motion will not significantly affect the determination of 
members. 

After membership selection a redetermination of the 
centre of mass could be made for the stars within, say, 
r < 10 pc (roughly corresponding to the tidal radius of 
the cluster) or r < 20 pc of the Hyades centre. Starting 
from the preliminary centre of mass only two iterations 
were needed to converge to corresponding determinations 
of the centre of mass and mass motion, and these are very 
robust in terms of the initial estimate. The results are 
given in the second and third lines of Table 3. We adopt 
a reference distance of 46.34 ± 0.27 pc, corresponding to a 
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Fig. 8. (a) Projected positions for the 218 candidate members, in Galactic coordinates (top); (b) projected velocity distributions, 
in Galactic coordinates (bottom). 
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Fig. 9. Projected velocities as a function of position for the 197 candidate members with available radial velocities. The residuals 
are given with respect to the velocity of the cluster centre in Galactic coordinates. HIP 13117 is omitted due to large errors in 
its space velocity. 



distance modulus m—M = 3.33±0.01 mag, for the objects 
contained in the Hipparcos Catalogue within 10 pc of the 
cluster centre (roughly corresponding to the tidal radius). 

The resulting value of the space velocity in equatorial 
coordinates (ICRS) is (-6.28, +45.19, +5.31) km s" 1 with 
a corresponding convergent point of (a, S) — (97?91, 6?66) 
for the inner 10 pc (134 stars), and (-6.32, +45.24, +5.30) 
km s _1 with a corresponding convergent point of (a, S) = 



(97?96,6?61) for the inner 20 pc (180 stars). The result- 
ing motion of the Hyades with respect to the LSR is de- 
rived from a solar motion of 16.5 km s _1 in the direc- 
tion (£,b) = (53°, 25°) (Binney & Trcmaine 1987), and 
is approximately (—32.7, —7.3, +5.9) km s _1 in Galactic 
coordinates. 

Although the results for the r < 10 pc and r < 20 pc 
samples are reasonably consistent, it should be evident 
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Table 3. Distance and velocity of the inferred centre of mass of the Hyades, as described in the text. The first line 
corresponds to the preliminary determination based on 142 stars in the central region. The last two lines correspond to 
the 'final' determination for the 134 stars within r = 10 pc of the cluster centre, and for the 180 stars within r = 20 pc 
of the cluster centre, respectively. 



Selection 


N 




be (pc) 






vc (km s 1 ) 




D 


V 






X 


V 


z 


u 


V 


w 


(PC) 


(km s" 1 ) 


Preliminary 


142 


-42.23±0.24 


+0.15±0.06 


-17.09±0.10 


-41.53±0.16 


-19.07±0.11 


-1.06±0.11 


45.56±0.27 


45.72±0.22 


r<10 pc 


134 


-43.08±0.25 


+0.33±0.06 


-17.09±0.11 


-41.70±0.16 


-19.23±0.11 


-1.08±0.11 


46.34±0.27 


45.93±0.23 


r<20 pc 


180 


-43.37±0.26 


+0.40±0.09 


-17.46±0.13 


-41.73±0.14 


-19.29±0.11 


-1.06±0.10 


46.75±0.31 


45.98±0.20 



that we are not in a position to provide an unambigu- 
ous value for the 'mean distance' of the Hyades, since the 
centre of mass is sensitive to the subset of stars used to 
calculate it, which in turn depends on the selection of stars 
contained in the Hipparcos Catalogue, as well as on the 
contribution of faint stars, white dwarfs, and secondary 
components of unresolved double systems. 

In principle the consistency of the results can be im- 
proved by using the redetermined centre of mass to carry 
out a new iteration of the membership selection. This pro- 
cess should ultimately lead to a consistent set of members 
and centre of mass. We used the centre of mass velocity 
for the stars in the inner 10 pc and a redetermination of 
the matrix in Eq. (17) from their space velocities, to carry 
out a second iteration of the membership selection. The 
result is that ten objects listed as (possible) members in 
Table 2 drop out as non-members. All these are located 
beyond 10 pc from the cluster centre listed in line 3 of 
Table 3. Hence, a redetermination of the centre of mass 
would lead to the same result for the stars in the inner 10 
pc and to convergence of the membership selection pro- 
cess. This illustrates the robustness of our membership 
selection procedure. 



jects with appropriate values of the z'E _1 z statistic but 
for which we have no radial velocity, and are thus unable 
to rule on the overall stellar space motion for that object. 
Future assignments can be made on the basis of new or 
improved radial velocities, or 7 velocities as new binary 
orbits are determined. In interpreting Table 2 it should be 
noted that objects with column x = but without known 
radial velocity are unlikely to be members, irrespective of 
their radial velocity. Objects with apparently 'reasonable' 
values of 7r and K- a d, but with large values of c (and hence 
column x = 0) will have discrepant proper motions (not 
evident from the table). 



6. Discussion of the Hipparcos parallaxes 

6.1. Parallaxes determined from the convergent point 

In Sect. 4 we demonstrated that the proper motions used 
by Schwan (1991) result in a different convergent point 
from that derived from the Hipparcos proper motion data, 
and that the Hipparcos proper motion data themselves do 
not lead to a unique convergent point, but to a successive 
selection of objects at each iteration of the convergent 



point method occupying a different location in velocity 
The positions of the resulting candidates in the Hertzsprung-J. e ( a closer study of the stars listed in Xable 2 an d our 



Russell diagram are shown in Fig. 7. Those candidates re 
taincd according to the additional criterion 7r > 10 mas are 
indicated separately. Fig. 8(a) shows the positions of the 
218 candidate members in Galactic coordinates. Fig. 8(b) 
shows the projected velocity distributions. Fig. 9 shows 
the projected velocities as a function of position for the 
197 candidate members with available radial velocity. We 
return to a discussion of these distributions in Sect. 7. 

In the final three columns of Table 2, we provide 
the distance, d (in pc), of each object with respect to 
the adopted cluster centre defined by 134 stars within 
r < 10 pc (see Table 3) (column v); the statistic z'£ _1 z, 
constructed with respect to the velocity of the preliminary 
centre of mass (column w); and the assignments '0', T', 
or '?' (column x) as our final membership indicator. A 
T' is assigned on the basis of the z'£ _1 z statistic alone 
(with the limits of c = 14.16, or c = 11.83 if the radial 
velocity is unknown, corresponding to the 3a confidence 
interval defined previously), and independent of distance 
from the cluster centre; '?' is assigned in column w to ob- 



final list of members reveals that no corresponding spatial 
bias is introduced). 

Before proceeding with a discussion of the resulting 
space and velocity distributions of the candidate members, 
we will now examine to what extent the distances inferred 
from the convergent point analysis using ground-based 
and Hipparcos proper motions respectively, are consistent 
with the Hipparcos trigonometric parallaxes. In principle, 
this should give additional confidence in the quality of the 
Hipparcos parallaxes and proper motions, and an insight 
into whether our explanations for the distance discrepan- 
cies resulting from previous ground-based proper motion 
investigations are correct. 

Fig. 10(a) shows the parallaxes inferred by Schwan 
(1991) from his proper motion analysis (and based on his 
published convergent point) compared with the Hipparcos 
trigonometric parallaxes for the same stars. Fig. 10(b) 
is again constructed using Schwan's proper motions, but 
using the space motion derived in Sect. 5.3 to calculate 
the parallaxes from his values of the proper motions. The 
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lower panels show corresponding results but now based 
entirely on the Hipparcos proper motions, using the con- 
vergent point determined by Schwan (c), and finally the 
space motion from Sect. 5.3 (d). The important feature 
of these diagrams is that Schwan's inferred parallaxes are 
systematically smaller than the Hipparcos trigonometric 
values, a trend which is visible in Fig. 10(b), and still not 
completely eliminated in Fig. 10(c), i.e. when the indi- 
vidual Hipparcos proper motions are used but taken in 
combination with Schwan's convergent point. In contrast, 
as seen in Fig. 10(d), there is no systematic difference be- 
tween the parallaxes derived from the Hipparcos proper 
motions and our present determination of the convergent 
point on the one hand, and the Hipparcos trigonometric 
parallaxes on the other. 



(a) 



ScbwFiTi s convergent point 



(o) 



(b) 



HIP convergent point 



(d) 



Fig. 10. This figure shows, for the stars used by Schwan (1991), 
the parallaxes inferred from the convergent point, the space 
motion and the proper motions, versus the the Hipparcos par- 
allax. The two left panels show the parallaxes inferred from 
Schwan's convergent point and space motion, using Schwan's 
proper motions (top) and the Hipparcos proper motions (bot- 
tom). The two right panels show the same, but now the space 
motion derived in Sect. 5.3 is used to calculate the parallaxes 
from the respective proper motions. 



80 



60 



^40 

CO 



20 



20 40 

7T (mas, HIP) 



60 



80 



Fig. 11 shows the parallaxes inferred from the Hipparcos| 

proper motions and the present determination of the space 
motion of the cluster centre of mass (ordinate) versus the 
Hipparcos parallaxes (abscissa) for all the kinematically 
selected members discussed in Sect. 5.4. For the subset 
of 7T > 10 mas objects, the excellent correlation between 
the two quantities indicates that the space motion and the 
trigonometric parallaxes are fully consistent. The outliers 
include stars with large errors in their astrometry, objects 
located close to the convergent point, or objects located at 
almost the same declination as the convergent point with 



Fig. 11. The parallax inferred from the Hipparcos proper mo- 
tions and the present determination of the space motion of 
the cluster centre of mass versus the Hipparcos parallaxes, for 
all the kinematically selected members discussed in Sect. 5.4. 
The outliers are stars with large errors in their astrometry, ob- 
jects located close to the convergent point, or objects located 
at almost the same declination as the convergent point with 
/i Q * 3> /i<5 (two other outliers fall off the top of the plot). The 
inset shows the relevant cluster region in more detail. 



Ha* ^> (which will lead to a greater probability of erro- 
neously accepting the object as a member). ^From Fig. 4 
the proper motions from Hipparcos appear systematically 
larger than those used by Schwan in both right ascension 
and declination. Applying a sign-test to the observed dif- 
ferences shows that this trend is statistically significant at 
the 95 per cent confidence level. 

We can reconcile these results as follows. The discrep- 
ancy between our reference distance modulus of m — M = 
3.33, and that derived by Schwan (1991) of m-M = 3.40, 
corresponds to 0.07 mag in distance modulus, or a ra- 
tio of 47.9/46.34 = 1.034 in distance. ^.From Eq. (2), 
ds/dn = |a»h|/|ms| x HW|Vh| x sin(A s )/sin(A H ), where 
H and S refer to parameters from Hipparcos and Schwan 
(1991) respectively. Using |V s |/|Vh| = 46.60/45.72 = 1.02, 
where Vh corresponds to the velocity of the 142 stars used 
for the membership selection (Table 3); and median values 
of sin(As)/ sin(An) = 1.003 and |/xh|/|ms| = 1-007 corre- 
sponding to the objects in common between the two de- 
terminations (and after transforming Schwan's data from 
B1950 to J2000) also leads to a combined difference of 
3 per cent. 

While the median differences between the two sets of 
proper motions are below 1 per cent, the derived clus- 
ter distance is also sensitive to the cluster's space veloc- 
ity - Schwan's larger radial velocity at the cluster centre 
(39.1 km s _1 compared to 38.6 km s _1 ), and larger angu- 
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lar distance from cluster centre to the convergent point, 
both lead to a larger space velocity, and to larger values 
of sin(A) for all stars. 

6.2. Comparison with previous parallax determinations 

In demonstrating that the Hipparcos parallaxes and proper 
motions together provide a consistent picture of the 
Hyades structure, space velocity, and dynamics, our re- 
sults provide independent evidence (in addition to that 
provided by the catalogue construction) that the trigono- 
metric parallaxes and their standard errors may be taken 
at face value. Since recent ground-based determinations 
of trigonometric parallaxes for candidate Hyades members 
have reached formal standard errors of a few milliarcsec, 
a comparison between these and the Hipparcos Catalogue 
values should therefore permit further insight into dis- 
crepancies between parallaxes (and resulting cluster dis- 
tance modulus) determined by different ground-based ob- 
servatories. In this section we undertake a first exami- 
nation of these differences. It is important to recall that 
the Hipparcos parallaxes are to be considered as absolute, 
while ground-based determinations require corrections to 
convert the measurements to absolute values, taking into 
account the parallax distribution of the reference stars. 

Table 4. Parallaxes for the objects in common between 
the Hipparcos Catalogue and Upgren et al. (1990). 



HIP 


van Altcna 


TTabs (HIP) 


Tabs (U 90) 


No. 


No. 


(mas) 


(mas) 


20485 


276 


21.08±2.69 


17.4±3.8 


20527 


294 


22.57±2.78 


18.5±4.8 


20563 


310 


19.35±1.79 


24.2±3.2 


20679 


363 


20.79±1.83 


23.4±4.4 


20827 


459 


17.29±2.23 


22.4±5.6 


20850 


472 


21.29±1.91 


29.5±5.2 


20978 


560 


24.71±1.27 


29.0±5.8 


21138 


645 


15.11±4.75 


22.0±5.0 



Eight of the 23 Hyades parallax stars observed by Up- 
gren et al. (1990), three of the stars from the list of Ianna 
et al. (1990), six of the 10 Hyades parallax stars observed 
by Patterson & Ianna (1991), the spectroscopic interfero- 
metric binary 51 Tauri observed by Gatewood et al. (1992) 
and Torres et al. (1997a), and the spectroscopic interfer- 
ometric binaries 70 Tauri and 78 Tau observed by Tor- 
res et al. (1997b, c), are also contained in the Hipparcos 
Catalogue. In addition, 60 of our candidate members are 
contained in the Fourth Edition of the General Catalogue 
of Trigonometric Parallaxes (van Altena et al. 1995). 

Mean magnitudes of ground-based parallax stars are 
typically about V = 11 mag, significantly fainter than 
the median of the Hipparcos programme, which leads to 
relatively large values of (J^hip ^2 — 4 mas for these 
comparison objects compared with the median value, of 
ov,hip ~ 1 mas, for the Hipparcos Catalogue as a whole. 

Fig. 12 shows the relationship between the Hipparcos 
parallaxes and the determinations derived from the Van 



Vleck observations by Upgren et al. (these are listed in 
Table 4, which also gives the correspondence between the 
HIP number and the van Altena number used by Upgren 
et al.). We have eliminated HIP 20605 (vA 334) from the 
comparisons because of its inaccurate Hipparcos parallax, 
largely as a result of its faint magnitude (V = 11.7 mag). 
Upgren et al. used the corrections from relative to absolute 
parallaxes according to the precepts of van Altena et al. 
(1991). The ground-based and Hipparcos parallaxes are 
in reasonable overall agreement, with (7Thip — 7Ttj9o) = 
— 2.5±4.6 mas, to be compared with the average correction 
from relative to absolute parallaxes, applied by Upgren et 
al., of 2.6 mas. 
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Fig. 12. The differences between the Hipparcos parallaxes, 
and 7T a bs determined by Upgren et al. (1990) from observa- 
tions made at the Van Vleck Observatory for the eight stars in 
common between the two programmes. 

Three Hyades objects observed by Ianna et al. (1990), 
and referred to there by their Johnson number (Johnson 
et al. 1962), were also observed by Hipparcos, and these 
yield 7r H ip = 19.17 ± 1.93 mas and 7ri 9 o = 18.9 ± 4.7 mas 
for HIP 20419 (HY 259); tthip = 21.10 ± 2.22 mas and 
7TI90 = 11.4±4.8 mas for HIP 21946 (HY 318); and tthip = 
18.44 ± 1.66 mas and tti 90 = 35.7 ± 5.3 mas for HIP 23498 
(HY 351, more commonly called vB187); the significant 
discrepancy for the latter may be partly attributable to 
the acceleration of the photocentric motion observed by 
Hipparcos. 

Fig. 13 shows the relationship between the Hipparcos 
parallaxes and the determinations derived from the Mc- 
Cormick and Stromlo observations by Patterson & Ianna, 
Tabs, gctp (these are listed in Table 5; note that their Ta- 
bles 1 and 3 refer to van Bueren numbers for the first four 
stars, and Johnson numbers for the last six stars). The 
comparison uses their reductions from relative to absolute 
parallaxes based on the statistical corrections determined 
from the General Catalogue of Trigonometric Parallaxes 
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Fig. 13. The differences between the Hipparcos parallaxes, and 
TTabs, gctp determined by Patterson & Ianna (1991) from ob- 
servations made at the McCormick and Stromlo Observatories 
for the six stars in common between the two programmes. To 
avoid overlap of the error bars, the Hipparcos values have been 
shifted by +0.2 mas for the McCormick comparisons, and by 
—0.2 mas for the Stromlo comparisons. 



(van Altena et al. 1991), with an average correction of 
2.8 mas. The corrections to absolute derived by Patter- 
son & Ianna (1991) on the basis of spectroscopic para- 
llaxes of field stars were slightly larger and result in an in- 
creasing discrepancy with the Hipparcos values for all but 
the McCormick observations of HIP 19834. On the basis 
of this small number of stars, the ground-based trigono- 
metric parallaxes (and especially those from the Stromlo 
observations) are seen to be reasonably consistent with 
those from Hipparcos, with reliable estimates of the stan- 
dard errors, and (7Thip — rrpigi) = —3.5 ± 3.5 mas for the 
GCTP-corrcctcd Stromlo values. 

Gatcwood et al. (1992) have used the Multichan- 
nel Astrometric Photometer at the Allegheny Obser- 
vatory to determine the parallax of the spectroscopic- 
interferomctric binary 51 Tauri (HIP 20087, vB24) of 
19.4 ±1.1 mas (employing a correction from relative to 
absolute parallaxes of 1.7 mas), compared with 7Thip = 
18.25 ± 0.82 mas. A more recent orbital parallax for the 
same object has been derived by Torres et al. (1997a), 
7i"orb = 17.9 ± 0.6 mas, a value which puts these two 



fundamental distance determinations of this object in ex- 
cellent agreement. Determinations for the binary 70 Tau 
(HIP 20661, vB57) by Hipparcos, tthip = 21.47±0.97 mas, 
and by Torres et al. (1997b), 7r orb = 21.44 ± 0.67 mas, as 
well as for the binary 78 Tau (HIP 20894, vB72, 9 2 Tau) by 
Hipparcos, 7thip = 21.89 ± 0.83 mas, and by Torres et al. 
(1997c), 7T rb = 21.22 ± 0.76 mas, are also in particularly 
good agreement, and may be taken as further evidence for 
the reliability of these separate distance determinations, 
as well as for the robustness of the Hipparcos parallax 
determinations for binary systems. 

It is noted that the subsequent determinations of 
the mean cluster distance by Torres et al. (1997a, b,c), 
based on these three specific objects, is sensitive to the 
proper motions adopted for the individual cluster mem- 
bers. While they derive a resulting distance modulus of 
3.40±0.07 based on 51 Tau, of 3.38±0.11 based on 70 Tau, 
and of 3.39 ± 0.08 based on 78 Tau, we can infer that any 
revised estimate of the mean cluster distance based on 
the use of the incrtially referenced Hipparcos proper mo- 
tions would yield essentially the same distance modulus 
as that presented here. For example, using the Hipparcos 
parallax for 51 Tau, our convergent point for the inner 
20 pc region, and the Hipparcos proper motions, wc find 
a mean cluster distance of 46.14 pc for the 53 stars used 
by Torres et al. (1997a), a value very close to our centre 
of mass value. The formula used by Torres et al. to de- 
rive distances (di = do(/j,o / in)(sm(Xi) / sin(Ao)), where the 
subscript zero refers to the reference object), utilises the 
PPM proper motion for 51 Tau (about 5 per cent larger 
than the Hipparcos value), and the PPM proper motions 
of the additional cluster stars, which are almost all smaller 
than those of Hipparcos (as discussed in Sect. 6.1). These 
effects together cause the systematically larger distances 
derived by Torres et al. for the cluster. 

McClure (1982) derived a dynamical parallax for the 
Hyades binary HD 27130 (vB 22, HIP 20019) leading to a 
distance modulus of 3.47 ± 0.05 mag, considerably larger 
than most astrometric determinations. The Hipparcos 
parallax, 7Thip = 21.40 ± 1.24 mas (distance modulus 
3.35 ± 0.12 mag) suggests that the distance modulus in- 
ferred by McClure was overestimated. The resulting M-L 
relationship for binaries (including vB 22) has been dis- 
cussed by Torres et al. (1997a). 
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As discussed by van Altena et al. (1993, 1994), and in 
the Introduction to the Fourth Edition of the General Cat- 
alogue of Trigonometric Stellar Parallaxes (van Altena et 
al. 1995), the heterogeneous nature of ground-based pa- 
rallaxes makes any comparisons between them and the 
corresponding Hipparcos parallaxes difficult to interpret 
in any unified manner. Thus, in establishing the system 
of the Fourth Edition of the GCTSP (van Altena et al. 
1995) three distinctly different problems were addressed: 
(1) the correction from relative parallax to absolute paral- 
lax; (2) the relative accuracy of parallaxes determined at 
different observatories; and (3) systematic differences, or 
zero-point differences between observatories. Whether the 
present results, and in particular Figs 12 and 13, suggest 
that the corrections from relative to absolute parallaxes, 
or the relative parallaxes themselves, have typically been 
slightly overestimated from ground-based observations, re- 
mains to be understood. 

A comparison for the 60 candidate Hyades members 
contained in the Fourth Edition of the GCTSP is shown in 
Fig. 14. Determination of the distance modulus based on 
104 Hyades members from the GCTSP yielded m — M = 
3.32 ± 0.06 mag (van Altena et al. 1997b), in good agree- 
ment with our present determination, suggesting that the 
GCTSP and Hipparcos systems arc indistinguishable to 
within the limits set by the ground-based parallax accu- 
racies, at least for the Hyades region. 

The are significant differences between the Hipparcos 
parallaxes (and proper motion components) and the in- 
dividual Hubble Space Telescope Fine Guidance Sensor 
observations reported by van Altena et al. (1997a). For 
the three brightest stars out of the four in common be- 
tween the two sets of observations (HIP 20563/vA310, 
HIP 20850/vA472, HIP 21123/vA627) the Hipparcos pa- 
rallaxes are between 26-42 per cent larger than the corre- 
sponding HST values. Discrepancies between the proper 
motion components reach 10-20 mas yr _1 . We offer no 
convincing explanation for these differences but, given the 
consistency of the Hipparcos measurements presented else- 
where in this paper, presently favour the Hipparcos values. 
Additional investigations will be required in order to sub- 
stantiate these claims. 

6.3. Lutz-Kelker corrections, and other effects 

Problems with the use of trigonometric parallaxes have 
long been recognised (Eddington 1913, Lutz & Kelker 
1973, Smith & Eichhorn 1996). One particular difficulty 
of interpretation arises because the distance estimate d = 
l/n is biased, with the error distribution of this estimate 
obtained by multiplying the probability density function 
of 7T by the Jacobian of the transformation from tt to d 
(Luri & Arenou 1997). As a consequence, with d = 1/tt 
used to estimate distances, an a posteriori correction is 
required - this statistical correction is applicable when 
determined for and applied to a particular sample popu- 
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Fig. 14. A comparison between the GCTSP parallax (van Al- 
tena et al. 1995) and the Hipparcos parallax (both in mas) for 
the 60 candidate members in common between the two cata- 
logues. 



lation. 'Lutz-Kelker' type luminosity corrections arise cor- 
respondingly, and may be understood qualitatively as aris- 
ing from the fact that the error volume beyond the dis- 
tance corresponding to the measured parallax is larger 
than the associated error volume at smaller distances; 
accounting for this effect makes all derived luminosities 
brighter by an amount depending on the relative paral- 
lax error, ov jix. Applicable corrections are not the same 
for a uniform space distribution of stars compared with a 
sample drawn from a concentrated population. 

The bias on the distance estimates for individual ob- 
jects in the Hyades cluster based on the Hipparcos para- 
llaxes can be calculated analytically (Smith & Eichhorn 
1996), depending on the relative error in the parallax, 
which is less than 0.1 for 78 per cent of the Hyades mem- 
bers (94 per cent have ct^/tt < 0.2). Hence for most of the 
members the expectation value of the bias in the distance 
will be less than 1 per cent. Since the same bias enters in 
the tangential velocities V Q * and V$, the process of select- 
ing members based on kinematics should not be affected. 
However, for individual objects the bias in v x and v z (de- 
rived from V a *,Vg and Vr) can be quite large (up to 20 
per cent overestimate for v x and 8 per cent underestimate 
for v z ). Thus particular care must be taken when inter- 
preting space velocities for individual objects with a large 
relative error on the parallax. 
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Smith & Eichhorn (1996) demonstrated that although 
one can calculate analytically the bias in the distance de- 
rived from the parallax of an individual star, the variance 
in the expectation value for the distance is infinite, so that 
a correction of the measured distance is not possible. How- 
ever, if the distribution of true parallaxes were known, the 
parallax of each individual object could be corrected, in a 
statistical sense, a posteriori. Given the observed parallax 
and the associated error, the most likely value of the true 
parallax, and hence the true distance, could be estimated. 
In principle one could deconvolve the distribution of mea- 
sured parallaxes (e.g. Lindegren 1995) to derive the error- 
free distribution of parallaxes. However, this necessarily 
assumes that the observed sample of stars is statistically 
representative of the underlying parent sample. In reality, 
the parallaxes are measured for a sample of stars subject 
to specific selection effects. 

In the case of Hipparcos the selection effects are com- 
plicated, and the parent sample of stars is known with 
only limited accuracy. In view of the uncertainties in the 
knowledge of both the parent sample of stars and the se- 
lection effects, and given that for the relevant values of 
OV/V < 0.1 individual magnitudes are likely to be biased 
at levels at or below about 0.01 mag (cf. Smith & Eich- 
horn, Fig. 6) no statistical corrections of the individual 
parallaxes have been attempted in the present work. Any 
resulting effect on the HR diagram, including the slope of 
the main sequence (the fainter stars, having larger errors, 
will have a larger bias) is expected to be small, and has 
not been considered further in this work. 

However, an assessment of the possible resulting bias 
in the mean Hyades distance and velocity, corresponding 
to a realistic space distribution of Hyades member stars, 
has been investigated by Monte Carlo simulation based 
on: (i) a synthetic cluster at the position of the Hyades, 
generated with an (albeit simplified) Plummer density dis- 
tribution with a core radius of 2.85 pc (see Sect. 8.1); (ii) 
the luminosity function for the Hyades derived by Reid 
(1993); (iii) a crude approximation to the Hipparcos se- 
lection of stars as a function of magnitude; (iv) errors on 
the parallax which increase with the apparent magnitude 
of the star. 

For 1000 realizations of the cluster, 220 stars were gen- 
erated and realistic errors reflecting the Hipparcos obser- 
vations (with a median increasing with apparent magni- 
tude) were added to the parallaxes. For each star the cor- 
responding (ICRS) values of b x , b y and b z were calculated, 
and the space velocity components were calculated from 
the true proper motions and radial velocities combined 
with the parallax. From these quantities the mean values 
of the components of b and v were calculated, as well as 
the distance and velocity of the cluster. 

These results suggest that our mean distance to the 
cluster and overall velocity are each overestimated by 
about 0.09 per cent. The coordinate transformations im- 
ply that the components of b are then overestimated by 



0.2, 0.06 and 0.03 per cent respectively, while the velocity 
component v x is overestimated by 0.8 per cent, v y by 0.06 
per cent and v z is underestimated by 0.29 per cent. These 
biases are typically small in comparison with the errors 
quoted for these quantities in Table 3. 

To these problems of working with trigonometric para- 
llaxes, we note finally one effect resulting from the specific 
observational configuration of Hipparcos, which results in 
a correlation length of 1-2° over which some small residual 
correlations will probably exist in the derived astrometric 
parameters (Lindegren 1989). Such an effect results in the 
error in the mean decreasing not as n 5 but probably more 
as n 35 , leading to a small underestimate in the final error 
on the estimated distance of the centre of mass. 

7. Structure and kinematics of the Hyades cluster 

7.1. Spatial distribution 

Fig. 8 shows clear evidence for a centrally concentrated 
group, possibly extended in b x (in the direction of the 
Galactic centre), with an evident correlation between the 
velocities in the Galactic coordinate directions y and z. 
Fig. 9 indicates velocity residuals which increase with dis- 
tance from the cluster centre, and which may be sugges- 
tive of systematic motions. In this section we compare 
these results with the expected space and velocity distri- 
butions predicted from iV-body simulations, and address 
the question of whether the cluster members display ve- 
locity residuals consistent with a co-moving system with 
constant space velocity, or whether there is evidence for 
systematic expansion (or contraction), rotation, or shear- 
ing motion due, for example, to the effects of the Galactic 
tidal field or passing interstellar clouds, or to the shear- 
ing effect of differential Galactic rotation. These effects, as 
well as the evaporation of stars through relaxation by stel- 
lar encounters, are considered to represent primary mech- 
anisms responsible for the disruption of open clusters. 

Realistic A^-body simulations of the dynamical evolu- 
tion of open clusters have been made by Terlevich (1987), 
de la Fuente Marcos (1995), Kroupa (1995) and others, 
aiming to reproduce features such as the observed distri- 
bution of cluster ages, mass density, binary distribution, 
and mass loss. Comparison with the models of Terlevich 
prove to be of particular interest since her A^-body inter- 
actions not only took into account specific forms for the 
initial mass function and mass loss due to stellar evolution, 
but were also supplemented by the effects of the Galactic 
tidal field and transient tidal shocks produced by passing 
interstellar clouds. Since it is known that the age distri- 
bution of Galactic open clusters barely extends beyond 
about 1 Gyr, and that disruptive encounters are likely to 
be responsible for cluster break-up, we may hope to in- 
vestigate whether the space and velocity distribution of 
objects beyond the tidal radius provides evidence for such 
a disruptive encounter. For example, Terlevich finds that 
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a cluster loses 90 per cent of its stars in 100 million years 
after a collision with a giant molecular cloud, in which the 
stars would acquire rather high velocities. 

Kroupa (1995) starts with a very large proportion of 
primordial binaries, and traces the stellar luminosity func- 
tion as a result of mass segregation, evaporation, and 
changing proportion of binary systems, assuming mod- 
els with (M) = 0.32 and M c = 128 M Q , near to the 
peak of the mass function of Galactic clusters. He takes 
the disintegration time to be the time taken for the num- 
ber density to reach 0.1 stars pc~ 3 , characteristic of the 
Galactic disc in the proximity of the Sun and, finding that 
these disintegration times are significantly longer than the 
lifetimes of real clusters (Battinelli & Capuzzo-Dolcetta 
1991), suggests that mechanisms other than internal dy- 
namical evolution must be responsible for cluster disin- 
tegration, such as impacts with giant molecular clouds. 
These encounters have been further modelled by Theuns 
(1992a, b), and confirm the systematic velocity signatures 
in the outer cluster regions, depending on the nature of the 
encounter, predicted previously (e.g. Binney & Tremaine 
1987, Sect. 7.2). 

The equipotential cluster surface becomes open, due 
to the effects of the Galactic tidal potential, at distances 
from the cluster centre referred to as the tidal radius r t 
(King 1962, Wielen 1974). These openings provide an es- 
cape route for stars to evaporate from the system. In the 
case of the Hyades, we find r t ~ 10 pc (see Sect. 8.3). The 
existence of an extended halo formed by stars outside this 
tidal radius, but sharing the proper motion characteristics 
of the rest of the cluster members was originally noted by 
Pels et al. (1975), and has since been reported for other 
clusters. 

A striking feature of Fig. 8 is that, although the tidal 
radius is of the order of 10 pc, about 45 stars are never- 
theless found between 10-20 pc, a result consistent with 
the simulations by Terlevich (1987) - these demonstrated 
that since the openings of the equipotential surface are on 
the cc-axis, stars can spend some considerable time within 
the cluster before they find the windows on the surface 
to escape through. Such iV-body simulation models con- 
sistently show a halo formed by 50-80 stars in the region 
between 1-2 tidal radii - some of these stars, despite hav- 
ing energies larger than that corresponding to the Jacobi 
limit, are still linked to the cluster after 300-400 Myr. 1 

These results appear to be weakly dependent on the 
slope of the IMF: Terlevich used a — 2.75, although de 
la Fuente Marcos (1995) has argued that for systems with 
large numbers of members the different models have very 
similar behaviour for the escape rate. [In any detailed in- 
terpretation of the spatial distribution of our candidate 



1 Another mechanism which can allow a cluster to hold 
on to stars that are beyond the tidal limit is binding by an 
angular-momentum-like non-classical integral, cf. Henon 
1970, Innanen et al. 1983. 



members in Fig. 8, it should be noted that certain se- 
lection effects operate in restricting the region of these 
diagrams which may be populated, due to the restricted 
region of a, (5 used for the study (Sect. 3) combined with 
the projection of the equatorial coordinates into Galactic 
coordinates.] 

It would seem natural to identify escaping stars, ev- 
ident both from numerical simulations and from our ob- 
servations, with the extended Hyades stellar group. The 
existence of such a system - stars having the same average 
motion as the cluster but a very much larger space and 
velocity distribution - was first suggested by Hertzsprung 
(1909), and followed up by Stromberg (1922, 1923), Eggen 
(1960, 1982 and references therein) and others. Recent 
models have been discussed by Casertano et al. (1993). 
This group is characterised by proper motions coincid- 
ing in direction with that of the Hyades, but considerably 
smaller in size. We will examine the relationship between 
this group and the central cluster in Sect. 7.2. 

A possible flattening of the Hyades cluster has long 
been debated. Flattening of the equipotential surfaces, 
perpendicular to and directed toward the Galactic plane, 
was predicted by Widen (1967) and in the iV-body mod- 
els of Aarseth (1973). Van Bucrcn (1952) noted that the 
cluster appeared to be flattened along the Galactic plane. 
Although this was not evident in the observational stud- 
ies of Pels et al. (1975), the suggestion re-emerged with 
a fairly clear indication of flattening in the outer region 
noted by Oort (1979), and subsequently by Schwan (1991). 
The effect is clearly evident in the A^-body simulations by 
Terlevich (1987, Fig. 8). 

A principal components analysis of the distribution of 
space positions for our Hyades members within 20 pc from 
the cluster centre shows that the cluster has a prolate 
shape. The major axis lies almost along b x in Galactic 
coordinates, making an angle of ~ 16° with the positive 
x-axis. The intermediate axis lies almost along b y , making 
an angle of ~ 16° with the positive y-axis. The short axis 
lies along b z . The axis ratios are 1.6 : 1.2 : 1, where these 
values have been derived from the standard deviations in 
position along the three axes, corrected for the median 
error in position along the corresponding axes (similar 
axis ratios of 1.5 : 1.2 : 1 are derived by multiplying the 
quartilcs of the distribution by the factor converting the 
quartile to the standard deviation for a normal distribu- 
tion). The inner 10 pc region of the Hyades is more nearly 
spherical. The fact that the shape of the outer parts of the 
cluster is prolate suggests that it is primarily extended in 
b x , although it is possibly also slightly compressed in b z , 
as shown in Sect. 8.2. This is consistent with the exten- 
sion being caused by stars slowly escaping through the 
Lagrangian points on the x-axis. 

Fig. 15 shows the cumulative distributions and num- 
ber densities versus distance from the adopted cluster cen- 
tre, for single stars and binaries for various mass groups. 
In the central 2 pc region, only stars more massive than 
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Fig. 15. In (a) cumulative distributions versus distance from 
the cluster centre are shown for single stars (solid line, 95 ob- 
jects); binaries (resolved or spectroscopic; dotted line, 75 ob- 
jects); and spectroscopic binaries only (dashed line, 58 objects). 
In (b) cumulative distributions are shown for M < 1M© (solid 
line 68 objects) ; 1 < M < 2M (dotted line, 97 objects) ; and 
M > 2M0 (dashed line, 26 objects). In (c) the corresponding 
number densities are shown for single stars (•, 89 objects) or 
binaries (o, 71 objects). In (d) corresponding number densities 
are shown for M < 1.2M© (•, 90 objects and M > 1.2M (o, 
91 objects). 



about 1 M Q are found (right-hand figures), and most of 
these are binaries (left-hand figures). This general effect 
was already noted by van Bueren (1952) and Pels et al. 
(1975), and is precisely as found in the simulations by Ter- 
levich (1987). We note in passing that the division between 
'hard' and 'soft' binaries, taken as the separation at which 
the (circular) orbital velocity (v = yjGM/r) is equal to 
the rms random motion is, for a system mass of 1 M 
and v — 0.25 km s _1 (see Sect. 7.2) is about 0.07 pc, or 
~ 5 arcmin at the distance of the Hyades. Numerical and 
statistical theories of star clusters predict the formation of 
a halo where the density falls as r -9 , with q between 3 and 
3.5 (King 1966, Spitzer 1975). Although for clusters with 
low central concentration such a power-law density distri- 
bution over the whole cluster is an inadequate represen- 
tation, from Fig. 15 we find q ~ 3.3 for single stars in the 
range r = 4 - 10 pc, q ~ 2.9 for M < 1.2M and q ~ 3.5 
for M > 1.2M Q . Thus less massive stars are more spa- 
tially extended than the more massive ones, with a lower 
density than for the high mass stars in the central regions, 
consistent with results from numerical simulations. Bina- 
ries follow the radial distribution according to their mass, 



implying that mass is the predominant segregation factor, 
rather than whether the system is binary or not. 

This mass segregation in turn implies that cluster lu- 
minosity functions derived from the central regions of the 
cluster will include an artificial flattening due to such mass 
segregation. The overabundance of bright stars is well 
established observationally, with Reid (1993) attributing 
this to mass segregation and stellar evaporation, while 
Eggen (1993) also shows that the luminosity function is 
depleted at the faint end compared with the observed field- 
star luminosity function. 

In the numerical models, massive stars sink to the 
central region, forming binaries, which become harder 
through the interaction with lighter stars, which in turn 
acquire enough energy to reach the outer parts of the clus- 
ter. There is a strong preference for energetic binaries to 
be formed among the heaviest members, which tend to 
segregate towards the dense central part of the cluster. 
Binaries themselves may not play a significant role in the 
evolution of the cluster, unless there is a significant popu- 
lation of primordial binaries (e.g. as concluded by Kroupa 
1995). Griffin et al. (1988) estimate that 30 per cent of the 
cluster members with 2.6 < My < 10.6 are radial velocity 
binaries. For systems brighter than My ~ 13 Eggen (1993) 
finds a photometric binary proportion of 0.4. Kroupa in- 
terprets his own simulation results as suggesting that the 
total proportion of systems in the central 2 pc sphere that 
are binary stars may be as high as 65 per cent. Such mod- 
els have been extrapolated by Kroupa (1995), Reid (1993), 
Weidemann et al. (1992) and others to derive a mass of 
around 1300 M Q (and around 3000 stars in total) at birth. 

Counting the stars with an SB or RV indication in col- 
umn s of Table 2 and those with an indication C, G, O, 
V, X, or S in column u as bona fide binaries, we find in 
our sample of members a binary fraction of 40 per cent. 
This fraction increases to 61 per cent for the stars located 
within 2 pc from the cluster centre, with almost all bina- 
ries in the central region being spectroscopic. 

7.2. Velocity distribution 

We now turn to an examination of the velocity distribu- 
tion within the cluster. As stated in Sect. 5.2 the internal 
velocity dispersion in the centre of the cluster is not re- 
solved with the accuracy of our present velocity data. The 
intrinsic dispersion expected for a cluster like the Hyades, 
in dynamical equilibrium, is ~ 0.2 km s^ 1 (van Bueren 
1952, Gunn et al. 1988), below the upper limit of the ob- 
served dispersion. 

Investigation of the possible systematic effects due to 
rotation of the cluster was investigated by Wayman et al. 
(1965), while Wayman (1967) estimated a limit on the 
contraction of K = —0.013 ± 0.015 km s _1 pc^ 1 . Hanson 
(1975) considered the possibility of expansion, contrac- 
tion, or rotation contributing significantly to the star's 
space velocities, and inferred that the convergent point 
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solution was sufficiently insensitive to rotation that ob- 
servable effects would be seen before impacting on the 
convergent point distance. Hanson also discounted signif- 
icant shear due to differential Galactic rotation. Gunn et 
al. devoted particular attention to the assessment of the 
velocity dispersion and the possible rotational flattening 
of the system. 

The residual velocities shown in Fig. 9 suggest that a 
number of stars in the outer regions of the Hyades show 
very substantial deviations from the mean cluster mo- 
tion. These deviations are, by definition, within '3cr' of our 
mean cluster motion, but a systematic pattern suggestive 
of a rotation or shearing motion in the cluster does seem 
to exist, although it is noted that the interpretation of 
a projection into two-dimensions of systematic structures 
in three dimensions is not necessarily clearly evident to 
the eye. Although other interpretations have been consi- 
dered (see below), the correlations between the residual 
velocities in Fig. 8 turn out to be fully consistent with the 
observational errors, as we shall demonstrate. 

In deriving space velocities for the cluster stars we 
make use of the observed vector (n, ^i a *,^S, Vr). This vec- 
tor is transformed to a space velocity, implicitly invoking 
a transformation to [V m ,Vs, Vr). On the assumption that 
the astrometric errors are uncorrelated the transformation 
of the observables to the vector (n, V ast , Vs, Vr) yields the 
covariance matrix: 

/ ° \ 

S ^ 

VO a%J 
With a = A v /n 2 , S is given by: 

-fi m acrl a 2 y?^o 2 + A v aa 2 af a 2 ^/j.scr 2 . 
-fisacr 2 . a 2 [i a ^[i & al a 2 fj,ja 2 + Ayda^ J 

Hence, even in the absence of correlations between astro- 
metric errors, the parallaxes and velocity components V m 
and Vs will in general be correlated. Moreover, because of 
the position of the convergent point of the system with 
respect to the cluster centre, \i a * is positive and [is is neg- 
ative for most cluster members, and hence the product 
Ha*Hs is negative. Thus for most stars the uncertainties 
in 7T and V m are anti-correlated, the uncertainties in n and 
Vs are correlated, and the uncertainties in V a * and Vs are 
anti-correlated, which will lead to systematic behaviour of 
the uncertainties in the sample as a whole. These system- 
atics will be transferred to the space velocities. 

To confirm the prediction of correlations in the ve- 
locity residuals we proceed as follows. If we assume that 
the motions of the cluster members are only due to the 
mean motion of the cluster, then the observations of the 
velocities of the cluster members will all have the same 
expectation value. We can then average all measured ve- 
locities to obtain a mean motion, with the uncertainty in 



the mean given by the mean of the covariance matrices 
of the individual members. This mean covariance matrix 
can then be used to construct the confidence region within 
which all residual velocities should lie. This is illustrated 
in Fig. 16 for all the Hyades members with a known radial 
velocity. The contours delineate the confidence region at 
confidence levels 68.3%, 95.4%, 99.73%, and 99.99%. From 
an eigenvector analysis of the mean covariance matrix, the 
minor axis of the distribution of velocity residuals is found 
to point in the direction I = 105°, b = 46°, which explains 
the flattened appearance of the distribution of residuals 
in Fig. 8b (rightmost diagram). 

An examination of Fig. 9 reveals a correlation between 
the velocity residuals (magnitude and direction) and the 
distances (parallaxes) of the stars. This effect, especially 
evident in the leftmost diagram, can be understood as fol- 
lows. The difference between the observed and true stellar 
parallaxes (An = 7r b s — TTtruc) is not correlated with the 
true parallaxes. However, adding An to 7Tt ruo implies that, 
on average, the stars with the largest observed parallaxes 
will have a positive An (and vice versa for the stars with 
the smallest observed parallaxes). So the sign of the par- 
allax error is correlated with the observed parallax. The 
correlation between An and V a * and Vs, discussed above, 
will then lead to a correlation between the observed dis- 
tances of the stars and the velocity residuals. Thus both 
the overall distribution of the velocity residuals, as well 
as the correlation of the direction of the residuals with 
spatial position, can be attributed to observational errors. 

Nevertheless, a large number of stars (32 out of 197) 
are located outside the '3cr' contour. About half of these lie 
beyond 10 pc from the centre of the cluster, with around 
one half of the 32 stars being binaries. This indicates that 
many of these stars may have suspect velocities. But even 
for the 165 stars inside the '3cx' contour, z'E X z is not dis- 
tributed according to the expected \ 2 distribution with 3 
degrees of freedom. A large fraction of the stars shows 
larger deviations from the mean motion, suggesting that 
the model 'mean motion plus error in the mean' is insuf- 
ficient to fully describe the residual velocities. 

To investigate the distribution of z'E _1 z further we 
restrict our attention to a 'high-precision subset' of the 
cluster members: stars without any indications of multi- 
plicity, with a radial velocity determined by Griffin et al. 
(1988), and with standard errors on the Hipparcos paral- 
lax and proper motions of less than 2 mas and 2 mas yr _1 
respectively. This selection results in a subset of 40 stars 
from which HIP 18962 is suppressed because it is located 
at a very large distance from the cluster centre (45 pc), 
and HIP 21788 is suppressed because of its large resid- 
ual v x (5 km s _1 ). Fig. 17 shows the cumulative distri- 
bution of z'E _1 z for these stars (lower solid line), where 
we have used the appropriate mean velocity of this subset 
v c = (-42.31, -19.08, -1.43) km s" 1 (in Galactic coordi- 
nates) to characterise the mean cluster motion. The dotted 
line shows the expected distribution in the absence of any 
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Fig. 16. Contour diagrams delineating the (projected) confidence regions corresponding to the mean covariance matrix of the 
space velocities for all candidate members with available radial velocities (cf. Fig. 8b). The confidence regions indicate the 
expected distributions of residual velocities in the absence of intrinsic dispersion. The contours correspond to the 68.3%, 95.4%, 
99.73%, and 99.99% confidence levels. The thick line, at 99.73%, is the c = 14.16 contour (see Eq. 16). The crosses correspond 
to the residuals between the observed space motions, and the mean (centre of mass) cluster motion. 



intrinsic dispersion in the velocities, from which it is clear 
that there is an extra dispersion unaccounted for in our 
model, or present in our observational errors. 

This extra dispersion may originate from underesti- 
mates of the standard errors in the radial velocities, from 
the Hipparcos astrometry, from the presence of undetected 
binaries (Mathieu 1985), from systematic motions (e.g. ro- 
tation), or simply from the intrinsic velocity dispersion 
of the cluster. If the cause of the extra dispersion is at- 
tributed solely to the quoted standard errors being under- 
estimates of the true external errors, then the standard 
errors in the astrometry would have to be increased by 
50 per cent (factor 1.5), the errors on the radial veloci- 
ties would have to be increased by a factor of 3, or both 
would have to be increased by 40 per cent. We consider 
such explanations unlikely. 

Given the fact that the observed binary fraction pro- 
vides a firm lower limit on the actual binary population, 
and that many workers consider that the true binary frac- 
tion may approach unity, the cumulative distribution of 
z'S 1 z has been constructed as before, but with a stan- 
dard error of 0.3 km s _1 added quadratically to each of 
the components v x ,v y , and v z (shown by the thick solid 
line in Fig. 17). The resulting agreement is now excellent. 
For a cluster mass of ~ 460 Mq , an estimated core radius 
of ~ 2.5-3 pc, and a half-mass radius of roughly 4-5 pc, 
the expected mean cluster velocity dispersion is 0.2-0.4 
km s _1 . Our 'high-precision subset' contains mainly stars 
in the central region of the cluster where the velocity dis- 
persion will be higher than the mean. Thus we postulate 
that the 0.3 km s _1 added to the standard errors of the 
space velocity components can be ascribed to a combina- 
tion of the internal motion of the cluster, possibly sup- 



plemented by the presence of undetected binaries in the 
high-precision sample which would contribute further to 
the observed dispersion. This is in good agreement with a 
recent estimate of the internal cluster velocity dispersion 
of 0.25 ± 0.04 km s _1 derived by Dravins et al. (1997) on 
the basis of a maximum-likelihood determination of the 
cluster's astrometric radial velocities. 

In principle the errors in the space velocities originat- 
ing only from the radial velocities can be decoupled from 
those in the (more homogeneous) proper motions by se- 
lecting objects in a thin 'parallax slice' in order to assess 
whether the cluster dispersion can be resolved exclusively 
in the proper motions. Restricting the data set to the par- 
allax slice 21-21.5 mas, and using only single stars within 
5 pc from the cluster centre, each proper motion is trans- 
formed into a component parallel to the direction of the 
convergent point (which can be corrected for the angular 
size of the cluster) , and a perpendicular component which 
should reflect only observational or intrinsic dispersions 
(the classical v and r components, e.g. Smart 1938). In 
the perpendicular direction the spread in the proper mo- 
tions as measured by the inter-quartile range is smaller 
than the median error. Although the spread in the paral- 
lel component is about twice as large as the median error, 
the parallel motion is very large (about 110 mas yr _1 ) 
so that even the 2 per cent spread in distance caused by 
the narrow parallax range selected leads to a considera- 
ble artificial spread in the proper motions, of the order of 
1-2 mas yr _1 . As a result, the intrinsic dispersion in the 
proper motions remains unresolved. 

Although our observations therefore appear fully con- 
sistent with uniform space motion, we also examined four 
possible causes of systematic structure in the space veloc- 
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Fig. 17. The cumulative distribution of z'£~ z for the high- 
precision subset of objects described in the text (lower solid 
line) compared to the expected distribution in the absence of 
intrinsic dispersion given by the \ distribution with 3 degrees 
of freedom (dashed line). The thick solid line has been con- 
structed by adding a standard deviation of 0.3 km s _1 to each 
of the components v x , v y , and v z . 

ities in an attempt to obtain a better fit to the residuals 
seen in Fig. 9 without resorting to the assumption of an 
unmodelled internal velocity dispersion: 

(i) systematic errors in the Hipparcos parallaxes or 
proper motions, or in the ground-based radial velocities, 
will affect the inferred velocity field. Geometrical consider- 
ations demonstrate that a spatially extended system par- 
ticipating in uniform space motion does not yield constant 
space velocities if the radial velocities, or fi a and/or 
are subjected to offsets of the type expected for our par- 
ticular observational quantities. For example, a constant 
offset in the radial velocity zero-point does not transform 
to a constant displacement in the space velocity because 
of the angular extension of the cluster on the sky. Simi- 
larly, a non-inertial 'spin' in the Hipparcos proper motion 
system (an effect whose elimination was the subject of 
considerable effort during the finalisation of the Hipparcos 
Catalogue) would lead to proper motion displacements of 
the form: 

A/j, a cos 6 = —uj x sin S cos a — ui y sin S sin a + ui z cos<5 
A/i5 = +uj x sin a — aj y cos a 

where (u) x ,u y ,u) z ) represents the non-inertial spin com- 
ponents of the proper motion system, and A/j, a cos 5 and 
A/15 are the resulting offsets in the individual proper mo- 
tions. We were not able to model both the magnitude and 
structure of the residuals as due to a combination of these 
zero-point errors, and large radial velocity offsets alone 



would lead to a very different structure in the residual 
velocities from that seen in Fig. 9. Taking only the mag- 
nitude of the residuals into account \w\ would have to be 
of the order of 10 mas yr" 1 , more than an order of magni- 
tude above the limits on the non-inertial spin-components 
of the Hipparcos proper motions noted in Sect. 3.1; 

(ii) objects beyond the cluster tidal radius will be 
subjected to systematic velocity perturbations from the 
Galactic tidal field, which will lead to a systematic pat- 
tern of residual velocities with increasing distance from 
the cluster centre, and periodic with time. The effects can 
be calculated on the basis of the epicyclic approximation 
(Binney & Tremaine 1987, Sect. 3.2). It is therefore ex- 
pected that the velocity residuals in the regions of the 
cluster beyond the tidal radius will deviate from a pattern 
of uniform space motion. However, the resulting velocity 
perturbations depend on the escape velocity and the time 
of escape, making it difficult to model them on the basis 
of only a few escaping members; 

(iii) the velocity residuals increase with distance from 
the centre, and superficially appear consistent with a gra- 
dient of Av ~ 0.3 km s _1 pc -1 out to distances of about 
10 pc, a gradient necessary to explain the largest velocity 
residuals in Fig. 9. Such a high value would exclude rota- 
tion as the source of the largest velocity deviations seen 
in Fig. 9, since this would require a mass orders of mag- 
nitude larger than the observed mass if the cluster is not 
to be disrupted on a short time scale. A lower rotation 
within the central 5 pc region is not required to explain 
the velocity residuals which, we re-iterate, are consistent 
with a non-rotating system and our observational errors; 

(iv) we have examined the possibility that the cluster 
recently experienced an encounter with a massive object 
causing a tidal shear in the outer regions of the cluster. 
The magnitude of space motion of the Hyades with re- 
spect to the LSR is rather large (34 km s -1 ), and since 
most known massive objects in the vicinity of the Hyades 
are molecular or atomic interstellar clouds having a much 
smaller motion with respect to the LSR, the consequent 
large relative velocity between the Hyades and any ob- 
ject encountered make the interacting system well suited 
to treatment using the impulsive approximation (e.g. Bin- 
ney & Tremaine 1987, Chapter 7, Theuns 1992a, b). In this 
approximation the relative velocity of the two colliding ob- 
jects is perpendicular to the velocity disturbances gener- 
ated in the less massive system. We can use the direction 
of the minimum velocity dispersion given by the results 
of the eigenvector analysis noted above to derive the di- 
rection of the relative velocity of a postulated encounter. 
Since the space motion of the Hyades is known, we can 
infer the minimum velocity with which the perturbing ob- 
ject moves with respect to the LSR. This turns out to be 
about 30 km s^ 1 , thus tending to exclude the possibility of 
a recent high-speed encounter with a nearby giant molec- 
ular cloud. Moreover, since the velocity increments im- 
parted to the perturbing object scale as 2GM/b 2 V, where 
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M is the mass of the perturbing object, V the relative im- 
pact velocity, and b the impact parameter, for any plau- 
sible values of b and V a very high mass, of the order 
of 10 6 M , is required for the perturbing cloud, uncom- 
fortably large compared with estimates for nearby giant 
molecular clouds (e.g., Dame et al. 1987). 

We conclude that systematic effects in our data, or ex- 
ternal perturbations of the velocity field, are not evident 
in our results, for which the parallaxes, proper motions, 
and radial velocities are all consistent with uniform space 
motion. In contrast, our results do not permit an unam- 
biguous or definitive assignment of cluster membership, at 
least beyond the cluster tidal radius, or as the velocity dis- 
crepancies between the individual 3-d space velocity and 
the mean cluster motion increase. 

Table 2 includes objects with low values of z'E _1 z, but 
with large distances (small parallaxes) for which the dis- 
tance from the cluster centre is large. Candidate 'escapers' 
can also be identified, on the basis of their small d, but 
discrepant motions. Of the list of 97 Hyades group stars 
listed by Eggen (1982) 92 are contained in the Hipparcos 
Catalogue and, of these, 8 are in the region of sky covered 
by the present study. Table 6 lists these objects, together 
with their distance from the cluster centre, the velocity 
deviation, and the value of z'E z. 

HIP 18692 and 26382 are among the 218 members al- 
ready contained in Table 2. HIP 22697 has a velocity con- 
sistent with membership, but with a parallax (9.27 mas) 
just below our adopted threshold, while HIP 18975 is only 
just outside our 3a membership limit. All of the stars 
in Table 6 are located beyond the tidal radius of 10 pc; 
HIP 18975 may be especially interesting in this context 
since it appears to be an example of an object emerging 
from the Hyades, still close to the tidal radius. HIP 21788 
(vB 110), rejected from our high-precision subset, may fall 
into the same category - while still within our 3cr kine- 
matical contour, it evidently has a somewhat discrepant 
motion, and was rejected as a plausible member on the 
basis of its radial velocity by Griffin et al. (1988). 

Breger (1968) used Stromgren photometry to classify 
false members of Eggen's Hyades-group list on the basis of 
their discrepant metallicity indices. Breger modelled the 
mctallicity distribution of the moving group as a 50:50 
mixture of Hyades stars with field stars of lower metal- 
licity. Breger's list contains five of the stars in our area 
of study, and two of these (HIP 13834 and 18692) are 
considered by us as members. The former does not oc- 
cur in the later Eggen lists, although it has the appropri- 
ate metallicity for membership according to Breger, and 
is well matched to the cluster kinematics despite its rela- 
tively large distance from the cluster centre (20.5 pc). The 
latter object also occurs in Table 6, is far from the cluster 
centre, and is only just contained within our kinematical 
selection limit. With a low metallicity according to Breger, 
further studies may well confirm it as a non-member. 



Table 6. Stars from Eggen's (1982) list of Hyades group 
stars, classified according to distance from our adopted 
cluster centre, d, Av, and z'E~ 1 z. 



HIP 


d 


Av 
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(pc) 


(km s ) 




12184 


21.2 


13.0 


22.09 


12189 


21.1 


15.8 


33.05 


12828 


25.7 


7.3 


18.25 


16813 


49.4 


19.4 


51.11 


18692 


45.9 


6.3 


13.67 


18975 


12.7 


5.5 


14.75 


22697 


63.4 


5.4 


3.32 


26382 


16.2 


3.0 


2.07 



While our present membership analysis is based only 
on kinematical arguments, mctallicity determinations clearly| 
contribute important additional information to the mem- 
bership, especially farther from the cluster centre. Mean- 
while, our findings underline a plausible connection be- 
tween the Hyades stars and escaping Hyades members. 

We summarise the various kinematical populations ev- 
ident from this study as consisting of the core, the corona 
(extending out to the tidal radius r t , Kholopov 1969), the 
halo (with r > r t but still dynamically bound to the clus- 
ter), and the moving group population (with r <; 2r t , and 
with similar kinematics signifying remnants of past mem- 
bership). 

8. Dynamical properties of the cluster 

8.1. Potential and density 

In the previous sections we have demonstrated confidence 
in the spatial distributions derived for the cluster mem- 
bers, and our objective now is to derive a description of 
the smoothed mass distribution within the cluster, in or- 
der to assess its dynamical behaviour, and its interaction 
with the Galactic potential. 

The form of the cluster potential determines the dis- 
tances at which it behaves dynamically as a spherical sys- 
tem or point mass. The potential at a position r is given 
by: 

N 



where the potential is expressed in units such that the 
gravitational constant is equal to 1, masses are expressed 
in solar masses, and distances are expressed in pc. Poten- 
tials were computed on a regular mesh of points within a 
20 pc radius sphere, centred at the centre of the cluster. 
To avoid local irregularities due to stars too close to one 
of the points where ^ was computed, the effect of a star 
was ignored if it was within 0.1 pc of that point (intro- 
ducing a softening parameter to model the potential at 
smaller separations has a negligible effect on the conclu- 
sions). The values on the mesh were interpolated to obtain 
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Fig. 18. The gravitational potential as a function of radial dis- 
tance, ^(r), for the objects within 20 pc of the cluster centre, 
constructed as described in the text 

the coordinates of points with specific rounded values of 
the potential. Finally, equipotential surfaces were deter- 
mined by a least-squares fitting of an ellipsoid through 
points with the same potential. 

Other than close to the centre of the cluster (r < 
1.5 pc), the ellipsoids could be approximated by spheres 
to better than 2 per cent indicating that, dynamically, the 
cluster has a high degree of spherical symmetry. Fig. 18 
shows the value of the potential as a function of distance 
r from the centre. For r > 9 pc, ^(r) is well represented 
numerically by: 

»(r) = 2 -f (21) 

where ^ is expressed in M Q pc^ 1 , assuming GM = 1 for 
the Sun. Beyond 9 pc from its centre the cluster can there- 
fore be approximated dynamically by a point mass, with 
a mass (corresponding to all of the Hipparcos stars consi- 
dered here) of mo = 240 M Q . This result is a consequence 
of the low space density of the stars farther out from the 
cluster centre. 

The results presented in Fig. 18 and Eq. (21) are re- 
presentative only of the 180 stars selected within the 20 pc 
radius sphere. If we make the assumption that the Hyades 
stars not contained in the Hipparcos Catalogue have the 
same distribution, and that the total mass within this 
sphere is M , the ordinates of Fig. 18, and the coefficient 
of Eq. (21), would have to multiplied by M /mo - the 
total mass of the cluster has been variously estimated at 
between 300 M (Pels et al. 1975, Oort 1979), 400 M Q 
(Gunn et al. 1988), and 460 M Q (Rcid 1992), according to 
the corrections made for duplicity, faint stars, and white 
dwarfs. 

To represent the mean density of the cluster, each 
star of mass m was replaced by a sphere with the same 
mass, but with a spatially extended and continuous den- 
sity distribution. The adopted model is of a constant den- 
sity within a sphere of radius s/2, linearly decreasing up 
to s, with s = 6 pc. Several density distributions were 
tested, all giving similar results provided that the radius 
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Fig. 19. (a) The mean density distribution within the cluster, 
o(r), in M pc -3 ; and (b) the mass within the sphere of radius 
r, in M Q , constructed as described in the text. 


"if the (sphere a is not significantly different from this value 



The same procedure as for the potential was then followed 
to determine equidensity surfaces. In practice, these were 
very close to spheres, confirming the spherical structure of 
the cluster out to a radial distance of about 7-8 pc. Be- 
yond this radius, the density is too small for the method to 
yield significant results. The resulting density distribution 
p(r) (in M Q pc~ 3 ) and the cumulative mass distribution 
(in Mq) are shown in Figs 19(a,b). As for the potential, 
the results scale as Mo/nio in order to take account of 
missing stars. ^From the values of the potential, the es- 
cape velocity of the star can in principle be determined 
although, as demonstrated in Sect. 7.2, in most cases the 
velocity errors are too large for individual conclusions to 
be drawn. 

The resulting density distribution has been compared 
to both a Plummcr model and a King model. For the 
former the best fit values of the core radius and central 
density are 2.9 pc and 1.8 M pc~ 3 , corresponding to a 
central velocity dispersion of 0.21 km s _1 . In the case of 
the King model, the best fit model has a core radius of 
2.6 pc, a central density of 1.8 M Q pc~ 3 and a value of 
*(0)/cr 2 = 2.6, where *(0) is the central potential (Bin- 
ney & Tremaine 1987, Sect. 4.4). This corresponds to a 
true central velocity dispersion of 0.24 km s _1 . 

For both models the central velocity dispersion is in 
good agreement with the value derived in Sect. 7.2. How- 
ever, results inferred from these particular models should 
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be viewed with some caution. The masses in both models 
fall short of the total mass observed in the Hyades. This 
is in part due to the fact that we are missing the faint 
end of the luminosity function in the sample under study. 
Adding fainter members of the Hyades may change the 
overall density distribution, while significant mass may 
also be present in the outer regions of the cluster. This 
is supported by the fact that the observed half-mass ra- 
dius (<~ 5.7 pc) is larger than the half mass radius for the 
models above (3.7 and 3.0 pc, respectively). 

Table 7. Moments of inertia (in M Q pc 2 ) with respect 
to the three coordinate planes (Galactic coordinates), r is 
the radius of the sphere (pc) and N the number of stars 
included. 
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20 
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8.2. Dynamical shape of the outer region 

The previous results concern primarily the central part of 
the cluster. The star density in the outer region is rather 
low, and these stars do not significantly modify the local 
potential, as demonstrated by the 1 jr form of the poten- 
tial in the outer regions. The dynamical approach to the 
study of the large-scale structure is through the compu- 
tation of the moments and products of inertia of the clus- 
ter. Table 7 gives the value of the moments of inertia with 
respect to the three principal coordinate planes, in Galac- 
tic coordinates, computed within spheres with increasing 
radii r. 

The evolution of the moments of inertia within 6 pc 
is rather chaotic, although demonstrating a tendency for 
a spherical dynamical shape. Beyond r = 8 pc, the mo- 
ments of inertia with respect to the y and z axes become 
prominent, confirming the extension of the star distribu- 
tion essentially along the x axis. The determination of the 
principal moments of inertia confirms this: the principal 
axis is very close to the x axis (within 10°), while stars lie 
both in the direction of the Galactic centre and anticentre. 

The moment of inertia with respect to the xy plane 
is systematically smaller than that with respect to the xz 
plane, indicating that the elongated part of the cluster is 
slightly flattened, and that the spread is slightly smaller 
in the direction perpendicular to the Galactic plane, as 
found by Oort (1979) and confirmed by Terlevich (1987). 

The angular momentum of the cluster, derived from 
the individual space velocities, is small and insignificant 
in the central region, but grows as more distant stars are 



taken into account, with a major axis tending towards the 
direction I — 125°, b = 50°. That the observed angular 
momentum is negligible for the main dynamical core of the 
cluster lends further support to the conclusions of Sect. 7, 
that the velocity residuals are a consequence of the obser- 
vational errors, and further studies would be required to 
demonstrate the existence or absence of significant angu- 
lar momentum or rotation. 

8. 3. The effect of the Galactic potential 

The cluster is immersed in the potential of the Galaxy, so 
that the equipotential surfaces, unlike the spherical sur- 
faces for isolated clusters, are distorted and eventually be- 
come open. 

In the disk, the location of these Lagrangian points in 
an open cluster can be calculated from the Oort constants: 



_ ( GM C \ 
X ^-\AA{A~B)j 



1/3 



(22) 



where M c is the total mass of the cluster (King 1962, 
Eq. 24) and A and B are Oort's constants. This distance 
is referred to as the tidal radius of the cluster, r t (al- 
though the volume defined by the equipotential surface 
is not spherical). Adopting A = 14.8 km s _1 kpc -1 and 
B = -12.4 km s" 1 kpc" 1 (Feast & Whitelock 1997) and 
using M c = 400 M Q gives r t = 10.3 pc. Since this dis- 
tance is beyond the limit where the cluster potential has 
the form of Eq. (21), we can conclude that the tidal radius 
is of order r t = 10 pc, the precise value depending on the 
value of M c . This is consistent with the consequences of 
the asymmetry in the evolution of the moments of inertia 
with r: beyond this distance the stars behave like com- 
panions of the cluster but under the predominant forces 
of the Galactic gravitational field. 

9. The Hyades Hertzsprung-Russell Diagram 

We have established well-defined spatial and velocity cri- 
teria for the assignment of cluster membership, at least 
within the central 10-20 pc region, and we now use 
these members to refine the Hertzsprung-Russell diagram 
with the objective of presenting a consistent picture of 
the observational (i.e. colour versus absolute magnitude) 
and theoretical (i.e. bolometric magnitude versus effec- 
tive temperature) relationships. Clusters provide an im- 
portant environment for testing associated stellar evolu- 
tionary theories, representing a number of stars which are 
considered, as a first approximation at least, to be at the 
same distance, co-eval, and of a constant metallicity. The 
Hyades, as the nearest moderately rich cluster, has been 
studied in detail for these reasons. 

While stellar models have been highly successful in 
matching the overall features of cluster colour-magnitude 
diagrams, unambiguous detailed model fitting has proved 
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more elusive. The major problems compounding these 
studies in the specific case of the Hyades have been the un- 
certain distance modulus of the cluster, and its associated 
depth, both conspiring to make the transformation to ab- 
solute magnitudes uncertain; the uncertainty in the mean 
metal content and in particular the discrepancies between 
photometric and spectroscopic determinations (e.g. Cayrel 
de Strobel 1982, 1990) leaving open the initial conditions 
for the stellar evolutionary models; the contribution to 
both coordinates of the observational HR diagram due to 
the contribution of (undetected) binary systems; and addi- 
tional complications associated with all such models such 
as the transformation from theoretical to observational 
quantities (requiring accurate bolometric corrections and 
colour conversions) and remaining theoretical uncertain- 
ties, primarily those associated with the theory of convec- 
tion, including the value of the mixing-length parameter 
and the possibility of significant convective overshooting 
(e.g. Maeder & Mermilliod 1980, VandenBerg & Bridges 
1984). 

Although the distance modulus of the cluster as a 
whole may have been assigned a small standard error in 
any given study, the depth of the cluster is such that the 
contribution to the intrinsic scatter of the observational 
main sequence from significantly different distances of the 
individual members may be substantial (see, e.g., Cayrel 
de Strobel 1982). The direct result of this depth effect is a 
spread in individual distance moduli of member stars lead- 
ing to a main-sequence population in the (V, B — V) plane 
less sharply defined than those of more distant clusters 
(cf. Figs 2 and 21). The consequence of this observational 
scatter is that it has been difficult to make a reliable es- 
timate of the helium abundance of the cluster members 
on the basis of model fitting which has, in turn, precluded 
unambiguous matching of evolutionary models to features 
such as the cluster turn-off. 

With the trigonometric parallaxes from Hipparcos, we 
are in a position to construct the observational HR di- 
agram with an accuracy of about 0.1 mag on individual 
values of My or Mboi (these errors are still dominated 
by the standard error of the parallaxes; in comparison, 
the error on V is typically 0.01 mag or smaller, and may 
be neglected). Our goals in this section are to derive an 
optimally constructed observational HR diagram for the 
cluster, and to use these observations in isochrone fitting 
to determine the cluster age. 

These goals are achieved in two successive steps: (i) 
using stellar evolutionary model fits to a carefully con- 
structed Mboi versus T e g diagram for the lower part of 
the main sequence, we will define the locus of the Hyades 
zero-age main sequence (ZAMS) and hence, through mod- 
elling, an estimate of the cluster's helium content; (ii) the 
observations over the complete part of the HR diagram 
in combination with theoretical isochrones will be used to 
determine the cluster age, using models with and without 
convective overshooting in the core. In both cases, it is 



important to suppress known binary systems from the ob- 
servational diagram, and to take into account the possible 
remaining biases in both coordinates caused by unrecog- 
nised duplicity. The influence of interstellar reddening ef- 
fects has been considered to be negligible for the Hyades 
stars (Crawford 1975, Taylor 1980) and is neglected in the 
following discussions. 

9.1. The Hyades ZAMS and He abundance 

The physical parameters (effective temperature, spectro- 
scopic gravity, and mctallicity) of a significant number 
of stars in the cluster have been derived over the last 
25 years, and with advances provided by recent solid-state 
detectors, excellent spectroscopic data are now available. 
We have selected 40 stars, with bolometric magnitudes in 
the range 3-6 mag, for which high- resolution, high S/N 
spectra are available (Cayrel de Strobel 1980, Branch et 
al. 1980, Cayrel et al. 1984, Cayrel et al. 1985, Boesgaard 
1989, Boesgaard & Friel 1990). From these studies, the 
metal content and effective temperature are known with 
high accuracy for each star, with T e ff determined to typi- 
cally 50 K, or even better for some objects. Further details 
of the determination of the [Fe/H] is given by Cayrel et 
al. (in preparation). 

The relevant data are listed in Table 8. Bolometric 
magnitudes were calculated from the V magnitude given 
in the Hipparcos Catalogue, the Hipparcos parallax, and 
applying the appropriate bolometric corrections of Bessel 
et al. (1997). The error in M v , and therefore Mboi, is 
still dominated by the error on the trigonometric paral- 
lax, rather than by the apparent magnitude (which can be 
seen from an inspection of the standard errors on the mean 
magnitudes for these objects given in the Hipparcos Cata- 
logue). Lutz-Kelker-type corrections have been ignored in 
view of the small values of a^j-K (see Sect. 6.3). Fig. 20 
shows the resulting positions in the (Mbd, logT e ff) dia- 
gram. For the subset of 20 stars for which no evidence 
of binarity is indicated in Table 2, error bars are given 
corresponding to the standard error on the trigonometric 
parallax contributing to the standard error in M^ a \. The 
errors on T e fj are harder to quantify. Although the less 
massive stars were all observed, reduced, and analysed 
with the same methods (observations and modelling) , the 
more massive stars were observed, reduced and analysed 
by different authors, with different model atmospheres and 
different effective temperature scales. We have assigned er- 
ror bars corresponding to ±50 K in T e g for logT c ff < 3.78, 
and to ±75 K for logT c ff > 3.78, the latter in part taking 
into account the effect of the higher rotational velocity. 

Table 8 yields a metallicity for the Hyades of [Fe/H] = 0.14±| 
0.05. The observational quantity, [Fe/H], the logarithm of 
the number abundances of iron to hydrogen relative to 
the solar value, is related to the metallicity Z, in mass 
fraction, through [Fc/H] = log(Z/X) - log(Z/X) Q for a 
solar mixture of heavy elements (X is the hydrogen abun- 
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Fig. 20. Hcrtzsprung-Russell diagram (Mboi, logT c ff) for the 
40 dwarfs listed in Table 8. The symbols follow the multiplicity 
flags given in Table 2: objects which are spectroscopic bina- 
ries or radial velocity variable are indicated '*'; objects which 
are resolved by Hipparcos or known to be double systems are 
shown as circles; one object (HIP 18658) with detected pho- 
tocentric acceleration, and one object (HIP 19504) possibly 
resolved in photometry, are shown by triangles. For the re- 
maining objects, error bars correspond to ±50 K in T e g for 
logTeff < 3.78 and to ±75 K for logT cff > 3.78, and to 
in Mboi- ZAMS loci as a function of mass were constructed 
as described in the text, and are given for the Hyades (dashed 
line) and solar (dotted line) metallicities given in Table 10. The 
location of the Sun is also shown. 



dance by mass). With the solar value (Z/X) Q = 0.0245 of 
Grevesse & Noels (1993a) and the adopted Hyades [Fe/H] 
we obtain Z/X = 0.034 ± 0.007 which is slightly, but 
significantly, higher than the solar value. This is the ob- 
servational quantity to be used in the models. The error 
of 0.007 on (Z/X) includes the error on the solar (Z/X) of 
about 11 per cent (Anders & Grevesse 1989). 



Earlier discrepancies in the determination of [Fe/H] 
for the Hyades, in particular the differences determined 
from photometric and spectroscopic observations, have 
been summarised by Cayrel de Strobel (1982, 1990). The 
present value is in reasonable agreement with (although is 
not independent from) more recent determinations, e.g. 
Cayrel et al. (1985) ([Fe/H] = ±0.12 ± 0.03, or 0.14- 
0.15 accounting for the difference in activity between the 
Sun and the Hyades dwarfs); VandenBerg & Poll (1989) 
([Fe/H] = ±0.15). Note that vB52 was rejected as repre- 
sentative of the cluster mean by Cayrel et al. (1985) on the 
basis of its colour anomaly, possibly related to the strong 
emission of Ha. 

The lower part of the Hyades main sequence relevant 
to this study is populated by low mass stars which are 
only slightly evolved, and therefore rather close to their 
'zero age' position, the zero-age main sequence (ZAMS) 
defining the locus on the HR diagram where the stars 
become fully supported by core hydrogen burning. Being 
of intermediate or low temperature, their spectra do not 
show He-lines, and their photospheric He abundance can- 
not be determined spectroscopically. Theoretical compu- 
tations of internal structure have shown that the ZAMS 
locus depends on the initial chemical composition, both 
in terms of mctallicity and He content. Comparison with 
theoretical models is therefore used to estimate the He 
abundance. 

In order to fit the data with improved theoretical cal- 
culations, new zero-age main sequence models have been 
calculated by one of us (YL) with the CESAM stellar evo- 
lutionary code (Morel 1993, 1997). Updated input physics, 
appropriate to the mass interval covered by this sample 
of Hyades stars (from about 0.8-1.6 M©) has been used: 
updated OPAL opacities (Iglesias & Rogers 1996) com- 
plemented at low temperatures by the Alexander & Fer- 
guson data (1994), nuclear reaction rates from Caughlan 
& Fowler (1988), a solar mixture of heavy elements from 
Grevesse & Noels (1993a) corresponding to the mixture 
used in opacity calculations and the CEFF equation of 
state (Christensen-Dalsgaard 1991). 

The mixing-length parameter used for all models was 
a = I /Hp = 1.64 pressure scale-heights (I is the mixing- 
length and Hp the pressure scale-height), derived from the 
calibration in radius of the solar model with the same in- 
put physics (see below). The same value was applied to 
the Hyades following investigations of visual binary sys- 
tems with known masses and metallicity (Fernandes et 
al., in preparation) resulting in similar values of a for a 
wide range of ages and metallicities. A value of a = 1.5 
was suggested by VandenBerg & Bridges (1984), although 
a somewhat higher value was also not ruled out. Models 
using standard mixing-length theory are very sensitive to 
the precise choice of a, although the region of sensitivity 
is restricted. For a summary of the effects of the major 
uncertainties associated with the use of such models, in- 
cluding the physical terms such as the reaction rates, the 
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Table 8. Data for the 40 stars for which high-resolution spectra provide accurate values for T cff and M bol for the 
metallicity determination of the main sequence. The columns are: (1) Hipparcos Catalogue number; (2) van Buercn 
number; (3) V magnitude from the Hipparcos Catalogue; (4) Hipparcos parallax; (5) parallax standard error; (6) 
absolute visual magnitude, M v ; (7) error on M v due only to the parallax error; (8) log of effective temperature; 
(9) bolometric correction from Bessel et al. (1997); (10) resulting absolute bolomctric magnitude; (11) [Fe/H]; (12) 
reference for [Fe/H]: BLT = Branch, Lambert & Tomkin (1980); B = Boesgaard (1989); BB = Boesgaard & Budge 
(1988); BF = Boesgaard & Fricl (1990); CCC = Cayrel, Cayrel de Strobel & Campbell (1985); CCS = Chaffee, Carbon 
& Strom (1971); F = Foy (1975). 
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opacities, the universality of the mixing-length parame- 
ter, and the relevance of overshooting, see Lebreton et 
al. (1995). In the sample of low mass stars considered 
here, stars are believed to be essentially homogeneous on 
the ZAMS so that modelling uncertainties associated with 
convective core overshooting are avoided. 

As the models strongly depend on chemical composi- 
tion (Y, Z in mass fraction) , a first set of calculations was 
made with a wide range of chemical compositions yield- 
ing a grid of theoretical ZAMS with different helium and 
metal contents. Interpolation between these ZAMS then 
yielded the value of the helium abundance giving the best 
fit to the low-mass stars for the mean observational value 
of (Z/X) given in Table 9. 



A final ZAMS (see Table 10) was computed for 
Y = 0.260 and Z = 0.024, corresponding to the mean 
Hyades metallicity. As shown in Fig. 20 this ZAMS fits 
the observational data rather satisfactorily, with the larger 
scatter for log Toff > 3.78 possibly originating from un- 
derestimated errors on the derived effective temperatures 
for the more massive stars, as noted above. We pro- 
vide for comparison, in Table 10 and Fig. 20, the cor- 
responding results which, with the same physical assump- 
tions, give a proper calibration of the observed lumi- 
nosity (L = 3.846(1 ± 0.005)10 33 erg s^ 1 ) and radius 
(R Q = 6.9599 10 10 cm) of the Sun at an age of 4.75 
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10 9 years 2 using Z = 0.0175 ± 0.0015 (Grevesse & Noels 
1993b). 

The resulting ZAMS for the Hyades lies significantly 
above that of the Sun. Although the values of Y found 
for the Sun (0.2659) and the Hyades (0.260) are rather 
close, the higher metallicity of the Hyades implies that 
the relative helium to metal enrichment ratio is smaller 
than that obtained for the Sun (for the relevance of the 
assessment of this enrichment ratio, see Pagel 1995). The 
value of Y = 0.26 is close to the value used by Vanden- 
Berg & Bridges (1984) and would appear to rule out the 
suspicions that the Hyades is helium deficient compared 
with field stars (Stromgrcn et al. 1982). A refined value of 
the He content of the Hyades will be presented in a forth- 
coming paper, combined with the analysis of a number of 
Hyades binaries, whose masses are very sensitive to the 
He content. 

The uncertainty on our final estimate of Y has two 
sources. The observational uncertainty on [Fe/H] gives 
an error on Y of 0.02, as shown in Table 9. Additional 
uncertainties arise from the determinations of T e g and 
Mb \. ^From three different ZAMS models computed with 
Z = 0.024, and with Y = 0.25, 0.26, and 0.27, we infer 
that an error of 0.1 mag on Mboi leads to an error of about 
0.025 on Y. For the lower part of the main sequence, and 
a given Mboi, an error of 50 K on T e g leads to an error 
of about 0.025 on Y. Assuming that our individual errors 
on Mboi and T e g are uncorrelated, we estimate that the 
observational scatter in the overall T e g j Mboi diagram re- 
sults in an error of only ±0.01 in Y. Combined with the 
independent error arising from the uncertainty in Z, we 
infer that the final error in Y is dominated by the mean 
value of [Fe/H] used for the Hyades, and is Y = 0.26±0.02. 

We note that the 'single' objects apparently located 
1 — 2a above the resulting main sequence have a metallic- 
ity close to the mean value of the 40 stars. They may be 
undetected binaries, although their velocity residuals with 
respect to the mean cluster space motion (see Fig. 9) show 
no evidence for systematic departures which might be ex- 
pected as a result of their perturbed radial velocity. Since 
these objects are not outliers in a colour-colour diagram, 
the T c ff values may also be slightly suspect. 

9.2. Isochrone fitting to the Hyades and the cluster age 

If the chemical composition of a cluster is known, the ob- 
servational HR diagram in combination with theoretical 
isochrones allow the cluster age to be determined. Fig. 21 
shows the HR diagram (My versus B — V) for 131 stars 



2 The age derived for the Sun by different authors ranges 
from 4.5 Gyr (Gucnther 1989) to 4.75 Gyr (Christenscn- 
Dalsgaard 1982). Guenther's comparisons of two solar 
models with ages of 4.5 and 4.7 Gyr did not yield signifi- 
cant differences in the derived solar helium content nor in 
the mixing-length parameter. 



Table 9. Correspondence between [Fe/H] and Y, Z and 
Z/X used for the theoretical models described in the text. 
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0.034 
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Table 10. log(L/L Q ) and logT c ff as a function of mass for 
the Hyades (Y = 0.260, Z = 0.0240) and solar (Y = 0.2659, 
Z = 0.0175) zero-age main sequences derived from the 
CESAM code, as described in the text (see also Fig. 20). 
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+0.018 


3.7615 


+0.138 


3.7832 


1.20 


+0.121 


3.7727 


+0.231 


3.7934 


1.30 


+0.301 


3.7932 


+0.400 


3.8124 


1.40 


+0.458 


3.8113 


+0.551 


3.8309 


1.50 


+0.598 


3.8286 


+0.687 


3.8534 


1.60 


+0.725 


3.8485 


+0.810 


3.8800 



where, in order to ensure minimal contamination from 
non-cluster members, objects have been retained only if 
they satisfy our kinematical membership criteria, and are 
drawn from the objects lying within the r < 10 pc radius 
of the cluster centre (Table 2). From this sample, we have 
eliminated a few with <tb-v > 0.05 mag. Error bars cor- 
respond to the standard error in My, estimated from the 
standard error in the parallax, and in B — V . The apparent 
magnitudes V, the B — V colours, and the standard errors 
in parallax and B — V were taken from the Hipparcos Cat- 
alogue. Among these stars 72, indicated as filled circles in 
Fig. 21, are not classified as (suspected) double or multi- 
ple (Table 2) nor variable (as compiled in the Hipparcos 
Catalogue). The open circles represent the double or mul- 
tiple systems, spectroscopic binaries, and variable stars, 
and these will not be used in further discussions of the 
main sequence modelling. The main sequence has signifi- 
cantly reduced scatter in comparison with the V, B — V 
diagram of Fig. 2, and may be compared with recent de- 
terminations from ground-based observations, e.g. Schwan 
(1991). 3 



3 Dravins et al. (1997) have been able to further reduce 
the scatter in the Hyades HR diagram by constraining the 
radial velocities and parallaxes of the members according 
to the hypothesis of uniform space motion of the cluster. 
Although a model-dependent approach, the reduced scat- 
ter in their HR diagram, especially towards the faint end, 
supports our conclusions about the absence of systematic 
velocity structure within the cluster, whilst confirming our 
estimate of the internal velocity dispersion, and providing 
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Table 11. Stellar rotation, vs'mi, for stars in the turnoff 
region. 



HIP 


M v 


B-V 


v sin i 








(km s ^) 


20542 


1.55 


0.15 


40 


20635 


0.85 


0.14 


75 


21029 


1.54 


0.17 


70 


21683 


1.23 


0.15 


115 


23497 


1.13 


0.16 


115 



As expected, most of the stars falling significantly 
above the main sequence are double or multiple stars. 
There is no sub-dwarf sequence, a fact already noted in 
previous work (Hanson & Vasilevskis 1983, Griffin et al. 
1988). One star, HIP 17962, is located below the main se- 
quence. This is V471 Tau (WD 0347+171), an eclipsing 
binary consisting of a K0V star and a white dwarf (Nelson 
& Young 1970). In the turnoff region (B — V < 0.2 mag), 
there are five stars with no indications of duplicity or vari- 
ability: HIP 20542 (vB47), 20635 (vB54), 21029 (vB82), 
21683 (vB108) and 23497 (vB129). The star with the 
bluest colour index is HIP 20648 (vB56), a known blue 
straggler (Abt 1985, Eggen 1995). 

After elimination of the identified binaries, the lower 
part of the main sequence has an increased scatter com- 
pared to that seen for the bluer stars, a fact partly at- 
tributable to the lower accuracy of the corresponding 
parallaxes (as illustrated by the error bars), but possi- 
bly partially related to chromospheric activity. Campbell 
(1984) has shown that many red Hyades dwarfs show 
colour anomalies which are found to correlate with var- 
ious indicators of chromospheric activity. The four third- 
magnitude red giants 7 (vB28, HIP 20205), S 1 (vB41, HIP 
20455), e (vB70, HIP 20889) and 6 1 (vB71, HIP 20885) 
Tau represent the cluster's giant branch. The stars vB41 
and vB71 are known spectroscopic binaries. The star vB28 
was reported to be double from the speckle result of Mor- 
gan et al. (1982), although this result was not confirmed in 
the speckle duplicity survey of the Hyades cluster of Ma- 
son et al. (1993), possibly because the star was observed 
when the seeing was approximately 2 arcsec. 

After accounting for known binaries and variable stars, 
the diagram still appears to be possibly contaminated 
by unrecognised binaries. In particular, the stars HIP 
20901 (vB74), 21670 (vB107) and 20614 (vB53) lie sig- 
nificantly above the cluster main sequence. HIP 20901 
and 21670 are apparently Am type stars (Abt & Morrcll 
1995), amongst which the high frequency of spectroscopic 
binaries is well established (Abt 1961, Jaschek & Jaschek 
1987); the possibility that these stars are binaries can- 
not therefore be ruled out. HIP 20614 is a fast rotator 
(vsini = 145 km s" 1 ) (Abt & Morrell 1995), and pho- 
tometry indicates a possible binary (Eggen 1992). These 



further evidence that the majority of binaries, at least 
with a not too large Am, have been identified. 



three stars have been omitted in the subsequent fitting of 
the main sequence. 

Rotation affects the colours of the stars, the effect de- 
pending on the equatorial velocity and on the inclination 
of the rotational axis with respect to the line of sight, lead- 
ing to shifts of a few hundredths in B — V and a few tenths 
in My, generally towards the red and to higher luminosi- 
ties. Values of vs'mi taken from Abt & Morrell (1995) for 
the stars in the turnoff region are given in Table 1 1 . 

In order to compare the observational HR diagram 
with theoretical isochroncs, evolutionary models were cal- 
culated using the same input physics as for the ZAMS 
models from which the initial helium abundance was es- 
timated (Y = 0.26, Z = 0.024, a = 1.64). Sequences 
were determined for masses of 0.8, 1.0, 1.2, 1.4, 1.7, 2.0, 
2.5, 3.0 and 4.0 M Q , from the ZAMS to the beginning 
of the red-giant branch. For each mass, two evolutionary 
sequences were calculated: a standard sequence, and one 
taking into account an overshooting of the mixed convec- 
tive core which extends the size of the convective core over 
a distance of 0.20 pressure scale-heights (the latter signif- 
icantly modifying the shape of the resulting isochrones 
and hence the estimated ages - the reference value of 0.2 
comes from Schaller et al. 1992). Representative evolution- 
ary states for each mass were extracted, and the Geneva 
isochrone program used to obtain isochrones with ages in 
the range 500-750 Myr in steps of 50 Myr, with and with- 
out overshooting. 

One of the main difficulties in the transformation 
of isochrones from the theoretical plane (Mb i, logT e ff) 
to the observational plane (My, B — V) lies in the 
colour/temperature transformation. As a temperature in- 
dicator, the B — V index has the disadvantage of also being 
sensitive to metallicity. Moreover, since the stellar surface 
gravity varies along the position on the isochrone, the in- 
fluence of gravity on the B — V index has to be taken 
into account. We adopted the calibration from Alonso et 
al. (1996) which allows derivation of the B — V colour as 
a function of T e g and [Fe/H]. This calibration, valid in 
the range 4000 K < T cS < 8000 K, was extrapolated to 
higher T c ff according to the results of Haywood (1997, pri- 
vate communication). The adopted calibration yields, for 
the Sun (T cff = 5780 K), a B - V index of 0.62, in good 
agreement with the recent estimation of 0.628 ± 0.009 by 
Taylor (1997). The transformed B — V colours were then 
corrected for the influence of gravity according to the re- 
lationships given by Arribas & Martinez Roger (1988). 
Finally, in order to estimate My from Mboi the V bolo- 
metric corrections from Bessel et al. (1997) were adopted. 

Fig. 22 shows the theoretical isochrones correspond- 
ing to 550, 600, 650, 700 and 750 Myr calculated with 
overshooting, superimposed on the observational HR dia- 
gram. Fig. 23 shows isochrones corresponding to 500, 550, 
600, 650, 700 Myr calculated without overshooting. On 
both figures the zero-age main sequence is also indicated. 
Our present results show that the five stars located in the 
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Fig. 21. Absolute magnitude, My, versus B — V , for the stars considered as reliable cluster members within r < 10 pc. Filled 



circles indicate objects which are not classified as (suspected) double or variable. Error bars correspond to the standard error 
in the Hipparcos parallaxes and B — V colour indices. Specific objects indicated are discussed in the text. 
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The 69 single stars with the location of the Z A MS 



and with isochrones corresponding to the range 550-750 Myr 
calculated with overshooting. The two objects indicated by dif- 
ferent symbols (arrowed) are discussed in the text. 



turn-off region, for which we have no evidence of duplic- 
ity, can be reasonably modelled either with an isochrone 
of 550 Myr without overshooting, or with an isochrone 
of about 650 Myr calculated with overshooting. As the 
turnoff region is sparsely populated, we have included in 
the HR diagram the star HIP 20894 (9 2 Tau, vB72). This 
is the brightest star in Fig. 21, not counting those in the 
giant branch, and is the brighter component of a wide vi- 
sual pair (with 6 1 Tau) and a well-known spectroscopic 
binary. A difference in V magnitude of 1.10 ± 0.01 and in 
B — V colour of 0.006 ± 0.005 between the components of 
the SB star were estimated by Peterson et al. (1993). We 
adopted these values in combination with the Hipparcos 
parallax to place these two stars in Figs 22-23, where they 
are indicated as different symbols (arrowed). 

Using this system as an additional constraint, it ap- 
pears that the primary component of the system (9 2 
Tau A) does not lie on the isochrones between 500 and 
600 Myr without overshooting, while an agreement with 
the isochrone grid including overshooting remains very 
satisfactory. Although not included in Fig. 22 to avoid 
crowding, we estimate that the theoretical isochrone cor- 
responding to 625 Myr, calculated with overshooting, pro- 
vides an optimum fit to the present observational data. 
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Fig. 23. The 69 single stars with the location of the ZAMS, 
and with isochrones corresponding to the range 500-700 Myr 
calculated without overshooting. The two objects indicated by 
different symbols (arrowed) are discussed in the text. 

For a given set of models (with or without overshoot- 
ing) we estimate an age uncertainty of about 30 Myr 
coming from the visual fitting of an isochrone to the 
observations. However, the cluster age determination is 
model dependent, and an uncertainty of about 15 per cent 
arises from the uncertainty on the amount of overshooting 
adopted in the calculations. Another significant source of 
uncertainty comes from the relationship adopted in the 
transformation between T e g and B — V, where the corre- 
sponding uncertainty in the age determination could reach 
some 20 per cent. The effect of rotation, relatively unim- 
portant in the Hyades (as shown in Table 11), is unlikely 
to lead to an overestimation of the age by 50 Myr (Maeder 
1971). From all of these considerations, and taking into ac- 
count uncertainties coming from the adopted models, the 
transformation between T e g and B — V, and the effect of 
undetected binaries, it is difficult to assign a very mean- 
ingful estimate of the uncertainty of our age determina- 
tion, which may reach 100 Myr, although an uncertainty 
of about 50 Myr may be a more realistic estimate. In sum- 
mary, our results suggest a cluster age of 625 ± 50 Myr, 
with observational evidence for the presence of convective 
overshooting. Support for overshooting from cluster main 
sequence fitting was already presented by Maeder & Mer- 
milliod (1980), who compared observational HR diagrams 
for 34 clusters in the age range spanning the Pleiades to 
the Hyades, and found that agreement could be obtained 
if the stellar convective cores are extended by a certain 
amount due to convective overshooting or other physical 
mechanisms. 

Previous estimates of the age of the Hyades from 
isochrone fitting ranged between 500 and 900 Myr (Barry 



et al. 1981). Cayrel de Strobel (1990) gave 655 Myr from 
the mean value of different age determinations. Recently, 
Torres et al. (1997a) found 600 Myr. Kroupa (1995) esti- 
mated a dynamical age of about 500 Myr. Direct compar- 
ison of our result with the isochrone-based ages quoted in 
the literature is complicated by differences in the models 
used. All of the previous estimates rely on models with so- 
lar composition, or interpolated from models having differ- 
ent metallicities. In this work, we have calculated models 
specifically for the Hyades abundance. 

We stress that distances to star clusters based on 
main sequence fitting to the Hyades must be corrected for 
chemical composition differences: thus the fitted m — M 
value for a system with solar abundance should be re- 
duced by roughly 0.13 mag to allow for the fact that 

[Fe/H] Hyades = 0.14. 



10. Conclusions 

The Hipparcos parallaxes and proper motions together 
provide a consistent picture of the Hyades distance, struc- 
ture and dynamics. They yield a cluster convergent point 
motion consistent with the individual trigonometric para- 
llaxes, and together explain the larger distance modulus 
derived from the most recent ground-based proper mo- 
tion investigations as originating from differences in the 
magnitude of the adopted cluster space motion, and small 
systematic effects in the ground-based proper motions. 
Conversely, the smaller distance modulus traditionally de- 
rived from a variety of ground-based trigonometric paral- 
lax programmes are attributed to errors in these ground- 
based parallaxes, a conclusion supported by the most re- 
cent distance modulus derived from consideration of the 
GCTSP parallaxes by van Altcna (1997b), in good agree- 
ment with our present results. There is good agreement 
with determinations using high-precision radial velocities 
(Stefanik & Latham 1985, Gunn et al. 1988). Recent dis- 
tance determinations to individual objects in the cluster, 
most notably the results of Torres et al. (1997a, b,c), are 
in excellent agreement with the Hipparcos trigonometric 
parallaxes, although their extrapolation to a correspond- 
ing mean cluster distance is again affected by systematic 
effects in the ground-based proper motions used. 

The combination of the Hipparcos astrometry with 
radial velocity measurements from ground-based pro- 
grammes provides three-dimensional velocities allowing 
candidate membership selection to be based on positional 
and kinematical criteria. A number of new cluster mem- 
bers have been found within 20 pc of the cluster cen- 
tre, and candidates can be classified as escaping mem- 
bers on the basis of their velocity residuals. No evidence 
for systematic internal velocity structure is found; rather, 
the results are fully consistent with a uniform cluster 
space motion with an internal velocity dispersion of about 
0.3 km s _1 . Spatial distribution, mass segregation, and 
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binary distributions are consistent with iV-body simula- 
tions. 

The cluster has a tidal radius of r t ~ 10 pc. Out- 
side this region, the stellar distribution is elongated along 
the direction of the Galactic centre and anti-centre, and 
is slightly flattened in the direction perpendicular to the 
Galactic plane. Inside this region, the cluster has spher- 
ical symmetry with a core radius of r c ~ 2.7 pc, and a 
half-mass radius of 5.7 pc. The presence of objects closely 
linked kinematically with the cluster core, but well beyond 
the tidal radius, probably originates from stellar encoun- 
ters and diffusion beyond the Lagrangian points. 

The well-defined observational main sequence has been 
transformed into a theoretical Mboi versus T c g diagram, 
from which fitting of the cluster zero-age main sequence 
yields a helium abundance of Y = 0.26 ±0.02. Theoretical 
isochrones matching the helium and metal content pro- 
vide observational evidence for convective overshooting, 
and yield a cluster age of 625 ± 50 Myr. 

With the caveat that the primary importance of the 
Hipparcos results is to provide individual distances to clus- 
ter members, rather than an estimated distance to the 
cluster centre of mass (a concept meaningful only in the 
restricted context of the cluster members contained in the 
Hipparcos Catalogue), our estimated distance to the ob- 
served centre of mass for the objects within 10 pc of the 
cluster centre is 46.34 ± 0.27 pc, corresponding to a dis- 
tance modulus m — M — 3.33 ± 0.01 mag. This mean dis- 
tance is, in practice, only marginally modified (formally by 
about 0.4 pc) for the derived centre of mass for Hipparcos 
objects within r < 20 pc of the cluster centre. 
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